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This  thesis  is  based  on  research  conducted  at  the  newly  constructed  Aeronautics 
and  Astronautics  Control  Sciences  Laboratory  at  the  University  of  Washington.  The 
set  up  of  the  UWAA  20-bay  planar  truss  is  described,  and  is  based  upon  a  similar 
experimental  test  bed  at  the  United  States  Air  Force  Academy.  Validation  of  existing 
open  loop  finite  element  models  of  the  truss  in  the  cantilevered  configuration  was 
conducted  for  this  research.  The  formulation  of  non-linear  and  approximate  linear 
models  for  the  two  primary  means  of  structuraJly-borne  actuation,  the  air  jet  thrusters 
and  the  reaction  mass  actuators,  is  also  explored.  Computer  non-linear  simulation 
and  experimental  testing  of  classical  analog  controllers  to  actively  suppress  vibration 
of  the  planar  bending  modes  of  the  truss  using  the  air  jet  thruster  and  reaction  mass 
actuators  independently  and  jointly  for  hybrid  linear/non-linear  actuation  is  reported 
in  this  thesis.  The  design  of  zm  optimal  controller  for  the  reaction  mass  actuators 
utilizing  linear  quadratic  Gaussian  techniques  is  reported,  and  includes  non-hnear 
simulation  of  hybrid  actuation  utilizing  the  optimal  reaction  mass  actuator  controllers 
and  the  classical  air  jet  thruster  controllers.  Recommendations  for  future  study  based 
upon  both  the  hcirdware  implementation  and  the  computer  simulation  are  suggested 
in  the  conclusions  section  of  this  thesis. 
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Chapter  1 

INTRODUCTION 


1.1  Background 

The  control  of  flexible  space  structures  is  of  paramount  importance  in  several  high 
visibility  space  programs  today,  including  the  Hubble  telescope,  Space  Station  Free¬ 
dom,  and  the  Strategic  Defense  Inititive.  All  of  these  systems  require  stable  platforms 
from  which  high  precision  optical  systems  may  be  employed.  Weight  restrictions  of 
establishing  orbit  dictate  the  use  of  lightweight,  rigid,  flexible  structures  as  basic 
building  blocks  for  these  space  platforms. 

Several  methods  have  been  proposed  to  deal  with  the  low  inherent  camping  in  large 
space  structures.  Passive  damping  through  the  use  of  viscoelastic  support  members 
in  the  lattice  work  of  the  structure  [6],  and  incorporation  of  shock  absorbing  isolators 
[22]  are  perhaps  the  most  efficient  answers  to  the  control  problems,  but  lack  the 
performance  of  active  control.  Increasing  the  mass  of  a  structure  helps  to  reduce 
the  amplitude  of  disturbances  and  lowers  the  frequency  of  structural  modes,  but 
runs  contrary  to  mission  requirements  of  light  weight  for  orbital  deployment.  Active 
control  of  beam  steering  mirrors  and  optical  lences  onboard  these  stations  [1,  17] 
sidesteps  the  problem  of  low  structural  damping  in  the  supporting  space  structure. 
However,  sometimes  low  frequency  structural  resonance  effects  may  be  uncontrollable 
by  the  high  frequency,  high  precision  mirror  and  optics  controllers.  Therefore,  we 
realize  that  some  degree  of  structural  stability  is  desirable. 

Space  realizable  actuation  (SRA)  is  an  additional  constraint  imposed  on  the  con¬ 
trol  of  space-based  systems.  Control  of  everything  from  airplanes  to  assembly  line 
robots  incorporates  one  simple  basic  premise  for  operation-there  is  a  sufficient  react¬ 
ing  medium  to  operate  against  for  control  purposes.  The  airplane  control  surfaces 
act  against  the  surrounding  air,  the  robot  is  securely  mounted  to  ground  in  an  as¬ 
sembly  line.  In  space,  the  luxury  of  a  reacting  medium  is  lost,  so  all  energy  created 
must  be  dispensed  with  by  onboard  systems.  SRA  is  of  two  main  forms.  The  first 
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is  a  non-renewable  type  realized  through  thrusters  and  rockets  capable  of  translating 
the  structure  in  space.  The  second  type  is  momentum  exchange  equipment  realized 
through  linear  reaction  mass  or  rotary  reaction  wheel  actuators.  While  not  capable 
of  translating  the  body  through  space,  these  types  of  SRA  may  be  electrically  driven 
and  thus,  using  power  supplies  such  as  solar  cells  or  nuclear  fission,  may  be  used 
without  worry  of  depleting  a  finite  quantity  of  fuel. 

Research  at  the  United  States  Air  Force  Academy  (USAFA)  [8]  points  towards  a 
unique  solution  incorporating  both  air  jet  thrusters  (AJTs)  and  structurally  borne 
reaction  mass  actuators  (RMAs).  The  AJTs  have  advantages  at  low  frequencies, 
high  amplitude  disturbances  typical  of  rigid  body  slewing  problems.  The  RMAs 
enjoy  a  higher  operating  bandwidth  than  the  AJTs  and  therefore  suggest  using  the 
two  methods  of  actuation  in  conjunction  with  one  another.  Hybrid  actuation  was 
the  terminology  coined  by  the  USAFA  team  when  researching  these  two  controllers 
working  in  unison. 

My  research  was  the  first  to  be  conducted  on  the  newly  constructed  planar  truss 
testbed  at  the  University  of  Washington,  Department  of  Aeronautics  and  Astronautics 
Control  Sciences  Laboratory,  which  duplicates  the  USAFA  research  planar  truss.  The 
goal  of  the  work  reported  in  this  thesis  was  to  re-examine  the  previous  work  on  the 
USAFA  truss  to  include  both  modeling  and  control  aspects,  establish  a  base-line 
controller  for  the  UWAA  planar  truss,  and  to  investigate  classical  multiple  input, 
multiple  output  and  optimal  design  techniques  in  hybrid  actuation. 

1.2  Problem  Definition 

This  thesis  reports  the  set  up,  modeling,  and  control  of  a  lightly  damped  planar  truss 
in  the  cantilever  configuration.  Two  methods  of  actuation  were  available  to  control 
the  structure:  the  air  jet  thrusters  and  the  structurally  borne  reaction  mass  actuators. 
Both  AJTs  and  RMAs  were  located  at  the  midpoint  and  the  tip  (free  end)  of  the 
truss  and  were  configured  to  control  lateral  motion.  The  location  of  the  actuators 
was  chosen  to  best  damp  the  second  bending  mode  of  the  truss,  and  duplicated  the 
original  USAFA  configuration.  The  accelerometers  were  also  located  at  the  tip  and 
midpoint  of  the  truss  and  sense  lateral  motions,  and  thus  provide  colocated  sensor- 
actuator  pairings. 


3 


The  goal  of  this  thesis  was  to  improve  upon  the  open  loop  stability  of  the  truss 
through  the  implementation  of  classical  and  optimal  controllers.  A  performance  goal 
of  10  %  modal  damping  is  suggested  by  Miller  and  Crawley  [15]  as  a  minimum  level 
of  damping  in  future  space  structures.  Modal  damping  and  time  responses  to  modal 
excitement  were  the  primary  means  of  measuring  success  in  attaining  this  goal. 

Measurement  of  all  frequencies  reported  in  this  thesis  is  in  radians  per  second 
(rad/sec).  Hertz  is  the  choice  of  structural  engineers  desiring  to  report  modal  fre¬ 
quencies,  and  has  merit  when  describing  the  response  of  structures  to  the  mechanical 
engineering  community.  However,  from  strictly  a  control  systems  design  standpoint, 
frequencies  reported  in  rad/sec  are  much  more  desireable.  For  instance,  the  conver¬ 
sion  of  frequency  responses  (Bode  plots)  to  easily  recognizable  transfer  functions  is 
easier  when  results  are  reported  in  rad /sec.  As  this  thesis  is  directed  towards  the 
controls  community,  the  rad/sec  frequency  convention  will  be  used  throughout. 

1.3  System  Configuration 

In  this  section,  the  hardware  and  analysis  equipment  used  in  the  UWAA  planar  truss 
control  experiment  is  described.  Computer  modeling  of  the  truss  and  actuators  will 
be  discussed  in  the  chapter  two.  An  appendix  is  also  included  which  describes  some 
of  the  idiosyncracies  and  special  considerations  of  the  equipment.  It  is  intended 
to  provide  aid  to  future  students  working  with  this  particular  hardware  but  is  not 
essential  to  the  academic  understanding  of  this  thesis. 

1.3.1  Truss 

Figure  1.1  show  the  experimental  set  up  of  the  truss.  The  planar  truss  is  23.2  feet 
long,  containing  20  square  bays  of  equal  size  with  alternating  diagonal  members  to 
increase  stiffness.  The  cantilevered  end  of  the  truss  is  bolted  to  a  rigid  steel  table 
which  is  in  turn  anchored  into  the  concrete  floor.  The  truss  rides  on  steel  ball  bearings 
to  minimize  the  friction  with  the  supporting  table.  The  lateral  members  of  the  truss 
are  affixed  to  rigid  steel  batans  to  increase  the  mass  and  hence  lower  the  natural 
frequencies  associated  with  the  lateral  bending  modes  [Ij. 
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Figure  1.2;  Truss  hardware  in  cantilever  configuration 
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The  first  four  bending  modes  of  the  truss  were  experimentally  determined  to  occur 
at  9.30,  55.8,  146.52,  and  262.32  rad/sec.  The  1/4  point  (station  5),  midpoint  (sta¬ 
tion  10),  3/4  point  (station  15),  and  truss  free-end  tip  (station  20)  were  the  locations 
kept  for  the  formation  of  finite  element  models  of  the  truss,  and  will  be  referred  to  in 
the  thesis  by  either  their  location  or  station  number  interchangeably. 

1.3.2  Air  Jet  Thrusters 

Two  pairs  of  air  jet  thrusters  are  mounted  on  the  truss,  the  first  pair  at  the  tip, 
the  second  at  the  midpoint.  The  pairs  work  together,  firing  transversely  in  opposite 
directions:  a  positive  voltage  from  a  control  signal  fires  one  of  the  pair,  a  negative 
signal  fires  its  opposite  twin.  A  regulated  high  pressure  supply  of  air  enters  a  manifold 
located  near  the  base  of  the  truss,  which  is  then  routed  up  to  the  four  AJTs  via 
separate  polyvinyl  tubes.  A  solenoid  valve  manufactured  by  MAC  valves,  model 
113B-551BA,  in  each  of  the  AJT’s  releases  air  from  the  source  tubes  through  the 
AJT’s  nozzles,  providing  thrust  which  may  be  used  alternately  to  excite  or  control 
the  transverse  bending  modes  of  the  truss. 

As  previously  noted,  the  AJTs’  locations  were  chosen  primarily  to  control  second 
mode  disturbances.  Other  locations  are  possible  and  warrant  investigation. 

1.3.3  Accelerometers 

Two  high-precision,  low-noise  Sundstrand  PA700  accelerometers  are  mounted  to  mea¬ 
sure  transverse  motion  at  the  tip  and  midpoint  of  the  truss.  Their  location  is  colocated 
with  the  actuators’  locations.  The  scale  factor  of  these  accelerometers  is  9.70  inches 
per  second  squared  per  volt.  The  accelerometers  have  a  useful  range  of  plus  or  minus 
15  volts.  Two  additional  accelerometers  are  mounted  on  the  reaction  masses  of  the 
RMAs. 

1.3.4  Reaction  Mass  Actuators 

Two  structurally-borne  reaction  mass  actuators  are  attached  to  the  truss  at  the  tip 
and  midpoint  lateral  degrees  of  freedom.  The  RMAs  are  Ling  Dynamics  V102  Shakers 
mounted  on  a  base  of  aluminum  construction  and  ride  on  ball  bearings  similar  to 
those  for  the  planar  truss.  The  shaker  DC  motor  is  connected  to  the  truss  node  via 
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a  stinger  which  allows  for  slight  misalignment  of  RMA  and  truss  node.  The  power 
supplies  for  the  RMAs  are  variable  gain  Ling  Dynamics  signal  amplifiers,  model 
PA-35.  An  accelerometer  is  rigidly  mounted  to  the  RMA  aluminum  housing  to  sense 
lateral  motion  as  noted  is  section  1.3.3,  and  provided  for  measurement  of  RMA  lateral 
motion  relative  to  the  truss  relative  motion  at  the  corresponding  connection  mode. 

A  photograph  which  shows  the  AJT  and  RMA  at  the  tip  of  the  planar  truss  is 
shown  in  figure  1.3. 


Figure  1.3:  Actuators  on  board  the  truss  tip 


1.3.5  Analog  Computers 

Two  types  of  analog  computers  were  used  to  implement  control  laws  in  this  experi¬ 
ment.  They  were  the  Electronic  Associates,  Inc.  EAI  TR-20  computer  built  in  1966 
and  the  more  recent  Comdyna  GP-6  analog  computer.  Despite  the  age  difference 
of  the  two  computers,  their  operation  and  performance  is  similar.  No  attempt  was 
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made  to  differentiate  between  the  two  computers  when  analog  control  circuitry  was 
implemented. 


1.3.6  Dynamic  Signal  Analyzer 

Data  acquisition  and  signal  analysis  was  performed  on  the  Hewlett  Packaxd  HP  35665A 
D3aiamic  Signal  Analyzer.  It  is  a  two-channel  FFT  spectrum/network  analyzer  with 
an  extended  frequency  range  well  beyond  our  requirements.  The  HP  35665 A  is  pri¬ 
marily  a  frequency  domain  analyzer,  but  was  also  used  for  time  domain  measure¬ 
ments.  The  analyzer  was  equipped  to  perform  both  random  noise  and  swept  sine 
device  characterization. 


1.3.7  Force  Gage 

Although  not  used  in  the  implementation  of  control  laws  in  this  thesis,  a  Transducer 
Techniques  force  gage,  model  MDB-10,  was  used  extensively  in  the  device  characteri¬ 
zation  of  the  AJTs  and  the  RMAs.  The  force  gage  has  a  linear  range  of  -10  to  10  lbs, 
and  is  calibrated  in  compression  only.  The  readings  in  tension  were  not  calibrated 
from  the  factory,  but  are  sufficient  for  device  characterization.  The  force  gage  TM-2 
amplifier  is  equipped  with  a  10  filter  capacitor  for  noise  suppresion,  resulting  in  a 
linear  response  for  signals  less  than  100  rad/sec.  The  frequency  response  for  the  force 
gage,  assuming  the  100  rad/sec  bandwidth,  is  calculated  and  shown  in  figure  1.4. 


1.3.8  Linear/Non-Linear  Analysis  and  Dynamic  Simulation  Package 

Computer  analysis  and  simulation  was  performed  using  SIMULAB/MATLAB  soft¬ 
ware  by  Mathworks,  Inc.  [19]  run  upon  DEC  workstations  equipped  with  X-windows. 
For  non-linear  simulation  the  ”linear”  numerical  option  was  used.  The  simulation 
package  was  designed  to  work  for  systems  containing  many  linear  blocks  with  a  few 
non-linear  blocks,  which  is  typical  of  the  hybrid  systems  simulated  for  this  thesis. 
The  algorithm  is  based  upon  the  work  presented  by  Frederick  and  Rimer  [5]. 
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Figure  1.4:  Force  gage  predicted  frequency  response 
1.4  Thesis  Outline 

Theoretical  development  and  experimental  verification  of  the  AJT,  the  RMA,  and  the 
UWAA  planar  truss  models  are  reported  in  chapter  2.  The  use  of  classical  control 
techniques  to  increase  eigenvalue  damping  is  discussed  in  chapter  3.  Both  linear 
model  evaluation  and  experimental  results  are  presented  in  chapter  3.  Chapter  4 
presents  an  optimal  controller  for  the  RMAs.  Conclusions  and  recommendations  for 
future  study  are  stated  in  chapter  5. 


Chapter  2 

SYSTEM  MODELING 


2.1  Overview 

Before  designing  a  controller  for  the  planar  truss,  a  viable  model  must  be  formulated. 
An  existing  model  for  the  planar  truss  in  the  cantilever  configuration  had  already 
been  developed  using  finite  element  methods  [8,  17]  ,and  requires  validation  for  this 
application.  Likewise,  the  controllers  had  been  previously  explored  [2,  8],  but  lacked 
the  experimental  implementation  that  forces  one  to  rethink  the  theoretical  models. 
This  chapter  re-explores  the  theoretical  development  of  the  models  used  in  this  the¬ 
sis  and  compares  them  with  the  actual  experimental  data.  Changes  are  made  to 
theoretical  models  to  better  reflect  the  experimental  data. 

2.2  Air  Jet  Thruster  Modeling 

Before  examining  the  truss  model  it  is  necessary  to  develop  models  for  the  AJTs, 
the  primary  means  of  excitation  and  control.  The  models  of  the  air  jet  thrusters 
may  be  mistakenly  overlooked  when  examining  the  frequency  characteristics  of  truss. 
This  error  may  be  extremely  serious  considering  the  inherent  non-linearities  and  the 
limited  bandwidth  of  the  AJTs. 

2.2.1  Theoretical  Development  of  the  AJTs 

The  AJTs  work  through  the  actuation  of  a  solenoid  valve  which  may  be  either  open  or 
closed.  Resulting  output  force  is  the  full-on  or  full-off  thrust  characteristic  of  AJTs, 
resulting  in  their  nickname  of  "bang-bang”  controllers.  The  AJTs  used  in  this  thesis 
are  connected  to  a  source  of  pressurized  air  at  65  psig,  and  each  develop  a  steady 
state  thrust  of  0.6  lbs  when  opened.  Ideal  on-off  control  is  not  possible,  however,  due 
to  the  dynamics  of  the  air  flow  in  the  jet  and  the  lags  of  the  solenoid  valve.  The 
relationship  between  the  output  thrust  f{t)  for  an  ideal  AJT  with  a  uniform  solenoid 
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valve  time  delay  tj  and  a  steady  state  output  force  to  a  control  signal  e{t)  is  given 
by  Hallauer  to  be  [8]: 

/(O  =  +  sgn{e{t  -  rj))]  (2.1) 

where 

s^n[e(<)]  =  1,  e{t)  >  ej 
sgn[e{t)\  =  0,  |e(0  <  cj 
sgn[e{t)\  —  —l,e{t)<ed 

and  c<f  is  a  small  ’dead  band’  voltage  which  prevents  the  AJTs  from  firing  for  small 
biases  and  noises  in  the  control  signal. 

Equation  2.1  is  valid  assuming  the  AJT  develops  full  thrust  immediately  after  the 
solenoid  opens  the  valve.  This  is  an  unrealistic  assumption,  however,  as  a  rise  time 
delay  will  occur  after  the  valve  is  opened  as  the  AJT  thrust  builds  to  its  steady  state 
value.  The  lag  is  caused  not  by  the  time  it  takes  for  the  valve  to  open,  but  by  the  time 
it  takes  for  the  force  output  of  the  AJT  to  rise  to  its  steady  state  value  once  the  valve 
opens.  The  valve  opening  may  therefore  be  considered  to  happen  instantaneously 
with  the  application  of  a  control  signal  commanding  it  to  open.  The  rise  time  of  the 
AJT  force  may  be  then  be  modeled  by  a  second  order  critically  damped  system  with 
natural  frequency  u;„.  The  natural  frequency  of  a  critically  damped  second  order 
system  is  equal  to  [16]: 


where  tr  is  the  time  it  takes  the  system  to  reach  95  percent  of  its  steady  state  value 
when  subjected  to  a  step  input.  By  modeling  an  AJT  ’back  to  back’  pair  acting 
in  opposite  directions,  and  writing  the  time  domain  step  response  of  a  second  order 
system,  the  following  equation  results: 

f{t)  =  F^ssgn[e{t)][l  -  exp(-u;„(t  -  5u;(t))(l  +  -  sw{t))]  (2.3) 

where  the  function  3w{t)  is  defined  as: 
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sw{t)  =  t,  when  e{t)  =  0  (2.4) 


time 

Figure  2.1:  AJT  signals 

Thus,  a  new  AJT  force  pulse  is  initiated  whenever  the  control  signal  passes  through 
a  zero  reference.  The  force  will  be  positive  when  the  control  signal  is  positive,  and 
negative  when  the  control  signal  is  negative.  A  picture  illustrating  the  operation  of 
an  AJT  with  a  control  signal  e{t),  a  switching  signal  s^n[e(^)],  and  AJT  force  output 
f{t)  is  presented  in  figure  2.1.  This  switching  is  realized  in  the  hardware  by  applying 
a  positive  voltage  to  the  transistor  of  the  AJT  to  be  opened.  The  transistor  then 
closes  the  power  signal  circuit  to  the  AJT  solenoid  which  opens  the  valve  and  allows 
pressurized  air  out  through  the  nozzle.  Likewise,  a  negative  voltage  command  signal 
to  the  AJT  transistor  will  open,  and  shut  off,  the  power  signal  to  the  AJT  solenoid 
and  allow  the  valve  to  close.  The  AJTs  are  set  up  to  work  in  opposing  pairs.  The 
opposite  AJT  operates  by  applying  its  force  when  a  negative  voltage  is  applied  to  its 
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switching  transistor,  and  shutting  off  the  force  with  a  positive  voltage. 

Equation  2.3  is  a  simple  equation  for  the  magnitude  of  the  output  force  of  the 
AJTs.  However,  due  to  the  non-linear  nature  of  the  AJTs,  it  is  important  to  note 
that  it  is  impossible  to  relate  the  force  of  the  output  signal  to  the  magnitude  of  the 
input  control  signal  e{t).  Theoretically,  a  sinusoidal  input  control  signal  of  0.1  volt 
magnitude  will  have  the  same  magnitude  output  force  as  a  10  volt  input  signal.  This 
limitation  is  caused  by  the  fact  that  the  transistors  in  the  AJT  circuitry  close  the 
power  circuit  for  all  voltages  above  a  small  threshold  value.  Despite  this  limitation, 
an  arbitrary  linear  input/output  relationship  is  assigned  for  the  design  problem.  This 
relationship  will  be  discussed  shortly. 

The  phase  of  the  AJT  is  fortunately  easier  to  describe  and  verify.  As  proposed  by 
Hallauer  [8],  the  phase  lag  <j)  of  the  AJT  due  to  the  time  delay  of  the  AJT  solenoid  is 
simply: 


<t>d  =  -U}Td 


(2.5) 


The  phase  lag  in  our  model  is  instead  better  described  by  a  second  order  system 
transfer  function.  Substituting  s  =  ju  into  the  second  order  transfer  function  [16]: 


Giju;)  = 


(s  -f  u;„)2 


s  =  ju; 


(jw  -f  u;„)2 


where  the  phase  angle  is  therefore: 


(2.6) 


<t>T  =  L 


(;u;-t-u;„)2 


(2.7) 


There  is  one  additional  factor  that  affects  the  phase  lag  of  the  AJTs  that  is  not 
the  result  of  the  actual  mechanical  operation  of  the  AJTs.  A  lead  filter  precedes  the 
transistors  in  the  solenoid  switching  circuit.  This  lead  filter  circuit  is  modeled  by  the 
transfer  function  [17]: 


G{s) 


2^ 


(2.8) 
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<l>  =  4>l  +  (f>r  (2.10) 

One  additional  point  should  be  made  here.  Because  some  linear  matrix  algebra 
packages  work  best  when  all  the  roots  of  a  system  are  ordered,  it  is  convenient  to 
depart  slightly  from  the  critically  damped  case  for  the  AJT  force  output  and  allow  a 
variable  damping  ratio  which  may  be  adjusted  to  be  close  to  one  (critically  damped) 
yet  just  under  to  allow  for  oscillatory  (unique-non-repeated)  roots.  In  order  to  do  so 
equation  2.6  may  take  the  following  form: 


G(s) 


s2  +  2Cu;„s  -I-  a;„2 


(2.11) 


where  it  is  recognized  the  the  damping  ratio  C  will  be  just  under  one.  The  results 
of  this  section  are  shown  in  a  frequency  response  in  figure  2.4.  In  section  2.2.2  an 
explanation  is  presented  as  to  how  the  scaling  on  the  magnitude  plot  was  determined. 
The  non-linear  block  diagram  for  this  system  suitable  for  simulation  is  shown  in 
figure  2.2; 


Figure  2.2;  AJT  non-linear  block  diagram 
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2.2.2  Experimental  Data 


Figure  2.3:  AJT  in  test  stand 

The  AJT  was  placed  into  a  test  jig  shown  in  figure  2.3.  The  configuration  allows  for 
the  measurement  of  AJT  output  thrust  when  the  solenoid  valves  are  opened  by  an 
electronic  control  signal  e(t).  The  rise  time  4  of  the  AJT  was  measured  by  comparing 
the  output  thrust  time  response  to  the  square  wave  power  signal  to  the  solenoids.  The 
rise  time  constant  t^  for  the  AJTs  averages  about  30  ms,  and  was  found  by  measuring 
the  length  of  time  the  AJT  required  to  produce  95  percent  of  the  steady  state  output 
force  after  the  solenoid  opened.  A  natural  frequency  for  the  second  order  rise  time 
was  then  found  to  be  =  200ms  using  equation  2.2. 
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Figure  2.4:  Normalized  magnitude  and  phase  plots  for  the  AJT 

The  AJT  frequency  response  was  obtained  by  exciting  the  input  transistor  switch¬ 
ing  circuit  with  a  random  noise  source  from  the  dynamic  signal  analyzer  from  the 
range  of  1  to  50  hz,  or  6.28  to  157.0  rad/sec  and  taking  the  Fast  Fourier  Transform 
of  the  output.  The  magnitude  portion  of  the  frequency  response  was  scaled  so  that  a 
one  volt  input  would  equal  a  0.6  Ibf  output.  This  is  an  arbitrary  scaling  as  the  the¬ 
oretical  development  (see  section  2.2.1)  shows  the  independence  of  the  force  output 
from  the  control  signal  input.  It  is  useful  to  assign  an  input/output  relationship  to 
the  magnitude  plot  in  order  to  make  a  linear  model  of  the  AJT  while  realizing  there 
will  be  error  in  the  linear  model.  Thus,  the  experimental  data  appears  to  have  less 
resolution  than  was  actually  attained.  Furthermore,  a  swept  sine  capability  was  later 
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added  to  the  dynamic  signal  analyzer  software  package  which  would  give  better  ex¬ 
perimental  frequency  information  than  the  random  noise  excitation  used  on  the  AJT. 
In  figure  2.4  the  experimental  data  is  represented  by  a  dashed  line,  the  theoretical 
linearized  model  is  shown  as  a  solid  line. 

The  shape  of  the  magnitude  plot  (figure  2.4)  confirms  the  validity  of  the  modeling 
performed  in  section  2.2.1.  It  is  important  to  note  that  the  second  order  behavior  of 
the  force  gage  filter  is  reflected  in  the  experimental  data,  and  is  likewise  added  to  the 
theoretical  linear  model  so  that  a  valid  comparison  may  be  made.  The  experimentally 
determined  response  closely  matches  the  theoretical  values  below  30  rad/sec  and 
remains  within  3  dB  at  56  rad/sec,  the  frequency  of  the  second  bending  mode  of 
the  truss.  The  force  of  the  AJTs  drops  off  significantly  in  the  regions  of  the  third 
and  fourth  bending  modes  of  the  truss  at  rates  of  40  db/decade  (figure  2.4  shows  a 
more  extreme  drop  off  rate  of  80  dB  per  decade  due  to  the  force  gage  filter  effects), 
illustrating  the  problem  of  using  AJTs  alone  to  control  the  higher  modes. 

The  phase  plot  of  the  AJTs  (figure  2.4)  is  free  of  the  problems  of  the  lack  of  a 
relationship  between  the  magnitude  of  the  input  signal  and  the  magnitude  of  the 
output  force.  The  AJT  phase  relationship  weis  clear  in  equation  2.10,  and  the  ex¬ 
perimental  data  confirms  the  phase  relationship.  Again,  it  is  important  to  note  that 
the  theoretical  linear  phase  information  has  been  run  through  a  second  order  filter 
representative  of  the  force  gage  to  match  the  experimental  data. 

2.2.3  Linear  Model  of  the  AJT 

The  AJTs  are  modeled  with  a  third  order  hnear  model  to  produce  the  plot  shown  in 
figure  2.4.  A  third  order  model  was  chosen  to  reflect  the  theoretical  development  in 
section  2.2.1.  The  third  order  model  maintains  the  lead  filter  in  the  AJT  circuit  as 
well  as  the  second  order  critically  damped  force  output  of  the  AJTs.  The  AJT  model 
is  given  the  arbitrary  magnitude  scaling  of  0.61bf/volt.  Defining  a  linear  model  for 
the  AJTs  does  not  allow  one  to  design  controllers  using  linear  techniques-the  results 
would  be  misleading  and  most  probably  not  realizable  in  the  hardware.  However, 
by  defining  a  linear  model  for  the  AJT  it  may  be  used  in  the  linear  model  packages 
for  assessing  the  performance  of  other  controllers  implemented  on  the  UWAA  planar 
truss.  Specifically,  vibration  suppression  may  be  measured  in  the  magnitude  portion 
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of  a  frequency  response  with  a  input  to  the  AJTs  and  output  measured  at  a  truss 
location.  Experimental  open  loop  data  may  be  scaled  against  the  open  loop  responses 
of  the  linear  model,  and  then  closed  loop  responses  may  be  compared  using  the  same 
scaling.  This  method  is  used  extensively  in  chapter  three.  Furthermore,  the  validation 
of  the  open  loop  model  may  be  performed  by  comparing  open  loop  linear  responses 
against  scaled  experimental  frequency  responses  as  is  done  in  section  2.5.3.  In  all  of 
this  testing,  phase  information  is  not  affected  by  the  nonlinearities  of  the  AJT  and 
may  be  compared  directly  between  experimental  results  and  linear  model  results.  The 
linear  model  transfer  function  is  given  by  the  AJT  equation 


40000s +  0.004 
^  381 .5s2  +  40570s  +  60000 


(2.12) 


2.3  Reaction  Mass  Actuators 

Another  means  of  space  realizable  actuation  for  vibration  suppression  is  implemented 
in  the  truss  control  problem  by  the  use  of  reaction  mass  actuators.  The  RMAs  in 
this  test  configuration  act  laterally  at  the  end  and  midpoint  of  the  truss,  although 
any  lateral  location  along  the  twenty  bays  of  the  truss  could  have  been  chosen.  Ac¬ 
celerometers  are  mounted  onboard  the  reaction  mass  in  order  to  measure  the  relative 
acceleration  between  the  truss  bay  and  the  reaction  mass  using  the  difference  between 
the  truss  bay  acceleration  signal  and  the  RMA  acceleration  signal. 


2.S.I  Theoretical  Development  of  the  RMA  model 

The  RMA  consists  of  a  moveable  reaction  (or  proof)  mass,  a  colocated  accelerometer 
sensor  on  board  the  proof  mass,  a  fixed  coil  to  which  a  current  is  applied  via  a  power 
amplifier,  and  a  magnet  inside  of  the  coil  to  which  the  structure  to  be  controlled  is 
attached,  and  a  supporting  flexible  suspension  system  which  suspends  the  magnet  in 
the  coil  and  connects  the  magnets  to  the  reaction  mass. 

The  theoretical  development  of  a  reaction  mass  actuator  model  was  investigated 
by  Zimmerman,  Horner,  and  Inman  [23].  However,  their  reaction  mass  actuators 
were  much  more  complex  than  the  type  implemented  here,  for  they  had  provisions  to 
feedback  velocity  and  position  signals  to  the  power  signal  to  the  coil  and  hence  create 
artificial  damping  and  stiffness  constants.  The  RMA  here  is  of  a  much  simpler  con- 
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struction  and  may  be  modeled  by  a  constant  coefficient  spring-mass-damper  system. 
A  schematic  of  the  RMA  is  shown  in  figure  2.5. 


The  force  produced  by  the  energized  coil  is  proportional  to  the  input  current: 

Mt)  =  BLiit)  (2.13) 

where  B  is  the  magnetic  flux  density  of  the  magnet  which  is  cut  by  the  coil  and  L 
is  the  length  of  the  wire  in  the  coil  as  a  function  of  coil  diameter  and  number  of  turns 
per  unit  of  coil  length.  The  product  of  B  and  L  is  specific  to  each  RMA  and  may  be 
defined  as  the  constant  Gi  with  the  units  of  Ibs/amp.  The  current  supplied  to  the 
coil  is  proportional  to  the  input  voltage  by  another  constant,  G^,  with  the  units  of 
volts/amp.  Using  equation  2.13  and  the  definitions  of  the  constants  Gi  and  G2  the 
force  produced  by  the  coil  is  given  by  [23]: 

Mt)  =  G,G2v{t)  (2.14) 

where  u(*)  is  the  total  voltage  applied  across  the  coil.  If  back  emf  damping  is 
neglected  and  there  is  no  equilibrium  position  restoring  voltage,  then  the  total  voltage 
across  the  coil  is  simply  the  command  voltage  e{t).  Equation  2.14  becomes: 
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m  =  G^G2e{t)  (2.15) 

The  assumption  that  back  emf  is  negligible  for  the  RMAs  in  this  paper  is  valid 
because  the  flexible  suspension  system  provides  for  much  greater  factor  of  damping 
than  back  emf  is  capable  of  generating.  However,  for  an  actuator  with  velocity 
feedback  on  the  proof  mass,  back  emf  would  be  essential  for  modeling  purposes. 


Figure  2.6:  RMA  simplified  system 

A  spring- mass-damper  system  representative  of  the  RMA  in  the  test  configuration 
(mounted  to  a  rigid  support  structure)  is  shown  in  figure  2.6.  The  damping  of  the 
flexible  support  structure  may  be  approximated  by  a  viscous  damping  constant  c. 
The  stiffness  of  the  support  structure  may  likewise  be  approximated  by  the  linear 
spring  constant  k.  Summing  the  forces  acting  on  the  reaction  mass  and  applying 
Newton’s  second  law  yields  the  following  relationship: 

F{t)  =  Fi,{t)  —  cx{t)  —  kx{t)  =  mx{t)  (2.16) 

ma$9 

Taking  the  Laplace  transform  of  equation  2.16  (assuming  no  initial  conditions) 
and  applying  the  results  of  equation  2.15  results  in: 


F{s)  =  GiG2E{s)  -  csX{s)  -  kX{s)  =  ms^X{s) 


(2.17) 
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Equation  2.17  may  be  rearranged  to  provide  a  relationship  between  input  voltage 
and  the  output  position  of  the  reaction  mass; 


X{s)  _  G1G2 
E{s)  ms^  +  cs  k 


(2.18) 


Since  F{s)  =  ms^X{s),  equation  2.18  may  be  rewritten  to  provide  a  transfer 
function  relating  output  force  to  input  command  voltage.  This  transfer  function  is 
important  to  test  and  calibrate  the  RMA  and  is  well  recognized  in  the  literature  as 
the  "stand  alone”  transfer  function  for  a  reaction  mass  actuator  [8,  23,  2].  A  typical 
frequency  response  is  shown  in  figure  2.7. 


Fjs)  ^  G,G,ms-^ 

E{s)  ms'^  +  cs  +  k 

When  the  RMAs  are  attached  to  the  truss  midpoint  and  endpoint  the  lateral 
dynamics  of  the  system  must  be  re-derived  to  include  the  dynamics  of  the  truss. 
This  was  performed  by  Barker  in  his  dissertation  on  control  structure  interaction  [2], 
and  is  repeated  here  from  first  principles.  A  diagram  of  the  RMA  on  a  truss  lateral 
attachment  node  is  shown  in  figure  2.8.  The  mass  of  the  RMA  magnet  is  assumed 
to  be  small  and  is  neglected.  For  simplicity  the  truss  is  modeled  as  a  spring-mass- 
damper  system  as  well,  realizing  that  from  this  form  it  is  simple  to  implement  into 
the  MAPMODES  finite  element  code.  The  MAPMODES  code  is  a  finite  element 
model  developed  by  Professor  W.L.  Hallauer  of  Virginia  Polytechnic  Institute  and  is 
formally  titled  Matrix  Algebra  Package/Structural  MODES.  It  is  used  to  model  the 
planar  truss  in  section  2.4. 

To  realize  the  equations  of  motion,  a  Lagrangian  approach  is  used  [14].  The 
variables  are  defined  in  figure  2.8: 


jC  =  T  -V  =  ^m^xl  +  ^m,xl  -  ^ka{xa  -  x,Y  -  ^k,xl 


(2.20) 


For  the  coordinate  x.; 


“  ka^^a  ^5)  kgXg 


=  m,x. 
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Figure  2.8:  RMA  aboard  truss  station 

d  dc 

I .  o  •  “  "niaXa 

dt  0X(i 

Qxan-c  ~  Ca{Xa  3?,)  ■  -  /fc(t) 

The  resulting  equations  of  motion  are: 


±dC_^ 

dtdqr  dqr  ~ 

TflgXf  ^a(^a  ^»)  d"  "1”  ^a(^'a  ®j)  “b 

"H  ^a{^a  ®j)  ~  C(,(x'o  Xj)  fbi^') 


(2.21) 

(2.22) 

(2.23) 


The  equations  of  motions  may  be  rewritten  as: 


m^Xj  “h  CjXj  ^a(Xa  Xg^  ^ai^a  ^3}  ”b  ^sXg  fbi,^^ 

nZoXa  "1“  Cd^Xa  Xj)  ~1“  ^a(X(j  X5)  ■  ^^(t) 


Which  lend  themselves  to  the  matrix  form  first  shown  by  Barker  [2]: 


(2.24) 

(2.25) 
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In  order  to  implement  the  above  matrix  equation  into  the  MAPMODES  finite 
element  program  it  is  necessary  to  first  designate  a  relative  coordinate  Xa,  =  Xa  —  x,. 
This  coordinate  transformation  removes  the  dynamic  coupling  in  equation  2.26  in 
favor  of  an  inertially  coupled  mass  matrix.  There  are  two  advantages  to  presenting 
the  system  in  this  manner.  One  is  the  ability  to  easily  incorporate  the  matrices  into 
the  finite  element  code.  Secondly,  the  eigenvalue  damping  appears  in  the  diagonal 
damping  matrix.  This  gives  us  the  hint  that  the  relative  velocity  between  the  reac¬ 
tion  mass  and  the  structure  will  be  an  important  element  of  our  control  strategies. 
Equation  2.26  written  in  terms  of  the  relative  coordinate  XaS  is: 


2.3.2  RMA  Experimental  Testing  and  Linear  Model  Approximation 

The  RMAs  to  be  connected  to  the  truss  tip  and  midpoint  were  individuallly  charac¬ 
terized  in  the  specially  fabricated  test  stand  shown  in  figure  2.9.  The  RMA  is  in  the 
’stand  alone’  configuration.  The  RMA  to  be  calibrated  is  attached  via  a  metal  stinger 
to  a  force  gage  which  is  affixed  to  a  rigid  wall,  and  rides  upon  the  ball  bearings  it 
would  normally  ride  upon  when  attached  to  the  structure.  The  ball  bearings  reside  in 
machined  grooves  to  preserve  RMA  alignment.  The  stinger  is  attached  to  the  RMA’s 
magnet  (which  is  suspended  in  the  coil  by  the  flexible  suspension  system),  connecting 
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Figure  2.9:  RMA  in  the  ’stand  alone’  experimental  test  configuration 

the  magnet  to  a  force  gage  designed  to  measure  axial  force.  The  stinger  is  equipped 
with  a  ball  socket  joint  at  each  end  to  insure  that  the  axial  alignment  of  the  RMA  on 
its  ball  bearings  is  not  affected  by  small  rotary  effects  of  the  truss  station.  Although 
not  necessary  on  the  test  stand,  the  stingers  are  considered  to  be  an  integral  part  of 
the  RMA  and  are  included  in  the  characterization.  Using  the  dynamic  signal  ana¬ 
lyzer,  a  swept  sine  characterization  was  performed  utilizing  the  signal  input  to  the 
RMA  amplifiers  while  the  force  gage  measured  the  output.  The  amplifier  was  set  to 
a  gain  of  2.  The  swept  sine  program  adjusted  the  signal  input  to  maintain  an  output 
force  level  of  approximately  0.5  Ibf  over  the  entire  frequency  range  of  the  test  without 
exceeding  2  Ibf  to  ensure  accurate  measurement  without  overstraining  the  force  gage. 

Figures  2.10  and  2.11  show  the  magnitude  and  phase  plots  of  the  tip  RMA  and 
midpoint  RMA,  respectively.  On  the  plots  are  the  experimental  frequency  responses 
as  well  as  the  theoretical  models  for  the  RMAs.  The  theoretical  models  are  augmented 


frequency  (rad/sec) 


frequency  (rad/sec) 

Figure  2.10;  Tip  RMA:  theoretical  (dashed),  experimental  (solid) 

with  a  second  order  filter  representative  of  the  force  gage  in  section  1.3.7.  Values  for 
the  mass,  damping,  and  stiffness  coefficients  in  the  theoretical  models  are  chosen  to 
reflect  the  experimental  data.  Close  agreement  between  the  RMA  models  and  the 
experimental  frequency  responses  is  achieved  below  100  rad/sec,  the  bandwidth  of 
the  force  gage  used  in  the  device  characterization.  This  validates  the  assumptions 
and  modeling  of  section  2.3.1  of  the  RMA  in  the  stand  alone  configuration.  The  final 
block  diagram  is  shown  in  figure  2.12.  The  experimentally  determined  coefficients 
are  listed  in  table  2.1. 

With  an  accurate  RMA  model  now  defined,  attention  may  be  turned  to  the  plant 
for  this  experiment,  the  20-bay  planar  truss. 


phase  (degrees)  Ibfi^volt  (dB) 


frequency  (rad/sec) 


Figure  2.11:  Midpoint  RMA:  theoretical  (dashed),  experimental  (solid) 


2.4  Truss  Modeling 


The  dynamics  of  truss  structures  are  well  understood  [14],  and  this  experiment  utilizes 
the  MAPMODES  finite  element  code  to  define  a  linear  state  space  model  for  the 
truss  developed  at  USAFA  [8]  and  modified  by  Barker  and  Schulthess  [2,  17],  In  his 
master’s  thesis,  Schulthess  reports  the  theoretical  and  experimental  validation  of  the 
finite  element  model  for  the  UWAA  planar  truss  used  in  this  controls  experiment. 
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F(s) 


Figure  2.12:  RMA  block  diagram 


Table  2.1:  Experimentally  determined  values  for  RMA  models 


constant 

RMA  tip 

RMA  midpoint 

G 

1  Ibf/volt 

1  Ibf/volt 

m 

0.0097255  lbf2^ 

0.0097542  Ibf^ 

tn 

c 

0.060  Ibf  sec/in 

0.065  Ibf  sec/in 

k 

16.5  Ibf/in 

18.0  Ibf/in 

41.12  rad/sec 

42.96  rad/sec 

C 

0.0750 

0.0776 

2.4-1  Truss  Theory 


The  planar  truss  may  be  modeled  as  a  simple  cantilevered  Euler-Bernoulli  beam  to 
provide  an  upper  limit  on  the  natural  frequencies  of  the  truss  and  the  approximate 
shapes  of  the  bending  modes.  Assuming  constant  bending  stiffness  El  and  constant 
maiss  per  unit  length  m  the  resulting  differential  equation  is  [14,  17]: 


d*y{x,t)  m  d^y{xj) 
dx*  ^  El  dt^ 


(2.28) 


This  equation  assumes  plane  sections  remain  planar,  thus  neglecting  shear  and 
rotary  effects,  as  is  the  case  for  the  planar  truss  testbed  where  analysis  is  limited  to 
only  the  first  four  lateral  bending  modes.  Limiting  the  truss  analysis  to  the  first  four 
bending  modes  is  not  arbitrary,  but  a  result  of  the  use  of  the  low-bandwidth  AJTs 
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to  excite  the  structure.  Assuming  a  separable  solution  and  applying  the  geometric 
and  natural  boundary  conditions  for  the  cantilevered  Euler-Bernoulli  beam  yields  the 
modeshape  equations: 


Yn{x)  =  An[{sin  k^yw^L  —  smh  ky/w^L){s\n  ky/iv^x  —  smh  kyjw^x) 

-t-(cos  ky/m^L  +  cosh  ky/w^L){cos  ky/m^x  —  cos  ky/w^x)] 

n  =  l,2,...,oo  (2.29) 

where 


k^ 


m 

El 


(2.30) 


and  Yni^)  is  the  modeshape  of  arbitrary  amplitude  for  a  beam  of  length  L  with 
natural  frequencies  defined  by  the  equation: 


cosh{w^kL)cos{w^kL)  =  — 1 


(2.31) 


2.4-2  Truss  Linear  Model 

A  low  order  linear  model  for  the  truss  was  obtained  by  using  the  MAPMODES 
program  and  the  model  reduction  technique  highlighted  by  Schulthess  in  his  master’s 
thesis  [17].  A  brief  outline  of  the  procedure  is  included  here  for  completeness. 

The  M.^PMODES  program  solves  the  generalized  63  degree  of  freedom  FEM  free 
vibration  problem 

[M\q  +  [K]q  =  0  (2.32) 

yielding  the  eigenvalues  of  the  full  order  system.  Since  only  the  first  four  bending 
modes  of  the  truss  systert  '-re  required  for  this  system  (as  they  are  the  only  lateral 
modes  of  the  system  within  the  frequency  range  of  the  primary  means  of  truss  dis¬ 
turbance,  the  AJTs),  modal  truncation  was  applied  and  a  damping  coefficient  was 
assigned  for  each  mode  (the  damping  coefficients’  experimentally  determined  values 
are  reported  in  section  2.5.2).  The  final  reduced  order  system  model  including  viscous 
damping  (assuming  no  cross-coupling  in  the  damping  matrix)  is  [17] 


[Mr\qr  +  {Cr]qr  +  [A'r]9r  =  fr 


(2.33) 
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Defining  a  state  space  model  is  possible  by  recognizing  that  equation  2.33  may  be 
rewritten  as; 

qr  =  -M-^Crqr  -  K^q,  +  MrVr  (2.34) 

The  state  space  model  is  therefore  of  the  form  (with  I4  the  4x4  identity  matrix  and 
O4  the  4x4  null  matrix): 


dt 


I4  O4 


(2.35) 


Where  qr  was  chosen,  in  the  case  of  the  model  with  the  four  retained  bending  modes, 
simlpy  the  displacement  of  the  tip,  3/4  point,  mid  point,  and  1/4  point,  such  that 
the  states  were  defined  as 


X  = 


/  \ 

Xi 

Qtip 

I2 

93/4 

9mid 

X4 

9i/4 

xs 

qtip 

Xg 

93/4 

Xr 

9mid 

<  > 

K  91/4  j 

(2.36) 


and  the  inputs  chosen  only  to  include  the  force  input  at  the  truss  tip  and  midpoint 
locations  such  that 


u  = 


«2 


(2.37) 


The  state  space  equation  then  has  the  form  of 


X  —  Asx  +  Bsu  (2.38) 

The  C  and  D  matrices  may  then  be  chosen  by  the  user  to  output  the  truss  accelera¬ 
tion,  velocity,  and  displacement  at  the  chosen  truss  locations.  The  numerical  values 
which  constitute  the  A  and  B  matrix  for  the  open  loop  truss  model  are  presented  in 
section  A.l. 
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Table  2.2;  Planar  truss  bending  mode  frequencies 


Modal  Frequencies  (rad/sec) 

Mode 

Beam  Theory 

FEM 

Truss  Hardware 

1 

9.60 

9.48 

9.30 

2 

60.44 

5.59 

55.80 

3 

169.14 

146.40 

146.52 

4 

331.19 

260.12 

262.32 

2.5  Open  Loop  Truss  Model  Validation 

Open  loop  testing  of  the  planar  truss  serves  several  purposes.  First,  it  shows  the 
agreement  between  the  both  theoretical  and  finite  element  models  to  the  existing 
hardware,  thus  verifying  the  modeling  techniques.  Second,  it  serves  as  a  basis  to 
fine  tune  the  linear  truss  model  to  reflect  the  discrepancies  between  the  theoretical 
calculations  and  hardware  test  results.  Finally,  there  is  no  analytical  method  for 
determing  the  modal  damping  of  the  truss  so  modal  damping  coefficients  had  to  be 
determined  experimentally.  The  truss  underwent  three  series  of  tests  used  to  confirm 
the  modeling  validity  and  help  with  the  tuning  of  the  model.  The  tests  were  open 
loop  acquisition  of  the  mode  shapes,  open  loop  modal  free  decay  time  responses,  and 
frequency  responses  of  mid  and  tip  accelerations  to  tip  excitation  via  the  AJTs. 

2.5.1  Mode  Shapes 

Mode  shapes  were  taken  by  driving  the  truss  at  its  first  three  bending  modes  using 
the  tip  mounted  AJTs  and  measuring  the  lateral  displacement  at  each  of  the  twenty 
truss  bays.  The  fourth  mode  created  displacements  below  the  threshold  of  the  test 
equipment,  consequently  it  was  not  measured  for  mode  shape  verification.  The  natu¬ 
ral  frequencies  for  the  truss  bending  modes*  are  shown  in  table  2.2.  The  mode  shapes 
were  normalized  to  the  point  of  greatest  displacement  and  graphed  along  with  the 

*  The  experimentally  determined  natural  frequencies  for  the  first  four  bending  modes  of  the  USAFA 
planar  truss  were  reported  to  be  9.80,  62.8,  155.2,  and  270.1  rad/sec,  slightly  higher  than  the 
experimental  values  shown  in  table  2.2  [8]. 
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theoretical  beam  theory  and  finite  element  model  predictions.  Stiffness  constants  in 
the  finite  element /linear  truss  model  were  changed  to  better  reflect  the  experimental 
data  natural  frequencies,  consequently  only  the  final  finite  element  model  results  are 
shown.  The  results  are  depicted  in  figure  2.13  for  the  first  four  modes  shapes  in 
relation  to  the  twenty  stations  of  the  planar  truss. 


Figure  2.13;  Normalized  bending  mode  shapes:  beam  theory  (dotted),  experimental 
(dashed),  and  FEM  (solid) 

2.5.2  Free  Decay  of  Bending  Modes 

After  the  mode  shapes  were  recorded,  the  next  step  was  to  determine  the  equivalent 
viscous  damping  coefficient  for  each  mode.  A  logarithmic  decrement  procedure  was 
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applied  to  obtain  the  viscous  damping  coefficients.  The  method  is  described  by 
Thompson  [20]  and  used  to  obtain  damping  coefficients  for  the  USAFA  truss  [18].  It 
is  outlined  here  for  completeness. 


A  viscously  damped  second  order  free  vibration  is  shown  in  figure  2.14  and  has  a 
a  solution  of  the  form: 

x{t)  =  Aexp"^‘^"‘sin{yi  -  Cu;„t  +  (f>)  (2.39) 

The  logarithmic  decrement  6  is  defined  as  the  natural  logarithm  of  the  ratio  of  any 
two  successive  amplitudes  [20].  It  is  formally  expressed  as: 


^  =  In  —  =  In  exp  sin(v/l  -  +  (f>) 

^2  sin(vT^^a;„(ti  +  ^)  +  <f>) 


(2.40) 


The  sine  terms  are  equal  for  multiples  of  — ,  therefore  equation  2.40  reduces  to: 


<5  = 


In 


exp 


exp  "d' 


(2.41) 
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Table  2.3;  Experimentally  determined  damping  ratios  ^ 


Mode 

(Cr) 

1 

0.0068 

2 

0.0042 

3 

0.0090 

4 

0.0160 

which  may  be  simplified  to  yield: 

^  =  In  exp^"”^  =  Cu;„—  (2.42) 

Given  that  the  damped  natural  frequency  is  equal  to  o^nv^l  —  0  [16],  equation  2.42 
becomes 

s/r=7 

The  damping  ratio  may  them  be  explicitly  defined  as 

c-  ,  ^ 

For  lightly  damped  structures  the  logarithmic  decrement  between  peaks  will  be  small 
and  the  expression 

i  =  l  (2.45) 

may  be  used  with  little  loss  of  accuracy.  A  similar  approach  for  peaks  separated  by 
n  cycles  will  show  that  the  logarithmic  decrement  may  be  expressed  as  [18]: 

6  =  iln—  (2.46) 

n  a:„+i 

The  logarithmic  decrement  for  each  mode  is  determined  experimentally  by  driv¬ 
ing  the  planar  truss  at  its  resonant  bending  frequencies  with  the  AJTs  until  a  full 
range  of  motion  is  recorded.  Then,  the  AJTs  are  shut  off  and  the  truss  is  allowed 
to  free  decay  until  no  motion  is  present.  For  the  first  mode  the  truss  endpoint  was 
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Figure  2.15:  Tip  free  decay:  actual  (top),  linear  model  (bottom) 


displaced  laterally  .5  inches  and  released  since  full  displacement  of  the  first  mode  by 
the  AJTs  could  cause  the  truss  to  shake  off  of  its  ball  bearing  supports.  The  accel¬ 
eration  at  the  truss  tip  and  midpoint  is  measured  during  the  decay  cycle.  Although 
the  damping  coefficient  (  was  derived  for  the  displacement  signal,  the  acceleration 
signal  will  decay  in  the  same  logarithmic  proportion  as  the  displacement  signal  for 
low  frequency,  lightly  damped  modal  vibrations  and  thus  is  sufficient  for  determining 
approximate  damping  ratios  for  the  truss  modes.  To  realize  an  average  logarithmic 
decrement  S  the  peaks  measured  are  at  90  percent  and  10  percent  of  full  displacement 
amplitude  and  formula  2.46  is  applied  where  n  cycles  are  counted  between  the  two 
measurement  points.  The  truss  midpoint  and  tip  values  are  averaged  as  well  for  all 
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Figure  2.16:  Tip  free  decay:  actual  (top),  linear  model  (bottom) 


modes  except  the  third  where  the  midpoint  lies  on  a  nodal  point  (see  figure  2.13). 
The  results  incorporated  into  the  model*  are  shown  in  table  2.3.  The  free  decay 

*The  damping  ratios  of  the  USAFA  truss  for  the  first  four  modes  were  determined  to  be  0.011, 
0.0022,  0.00278,  and  0.00632  [18].  The  modal  damping  of  the  USAFA  truss  was  lower  for  all 
modes  except  the  first  mode,  which  is  about  60%  higher  for  the  USAFA  truss  compared  to  the 
UWAA  truss.  The  damping  ratios  for  the  UWAA  truss  were  calculated  before  the  final  test 
instrumentation  and  hardware  was  placed  upon  the  truss  so  that  theoretical  calculations  and 
design  could  be  performed.  Later  testing  of  the  truss  indicates  that  a  value  of  damping  in  the  first 
mode  is  closer  to  0.012,  which  closely  reflects  the  Academy  data.  Although  0.0068  is  the  value 
used  throughout  this  paper,  special  mention  will  be  made  to  compare  the  damping  to  the  more 
recent  data  acquistion.  The  experimental  data  in  figure  2.15  has  been  changed  to  reflect  the  more 
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plots  of  the  experiment  data  and  the  linear  FEM  are  shown  in  figures  2.15  and  2.16. 
It  is  important  to  note  that  this  development  assumes  only  linear  viscous  damping, 
and  therefore  is  subject  to  errors  caused  by  non-ljnear  Coulomb  (dry-friction)  damp¬ 
ing  and  structural  (hysteretic)  damping  which  comprise  the  main  energy  dissipation 
features  of  the  planar  truss. 

2.5.3  Open  loop  frequency  response 

The  open  loop  frequency  response  of  the  truss  model  is  obtained  using  the  tip  mounted 
AJTs  to  excite  the  structure.  The  swept  sine  dynamic  signal  analyzer  package  is  used 
with  an  input  voltage  of  plus  and  minus  one  volt.  The  acceleration  signal  at  the  tip 
of  the  truss  is  measured.  A  frequency  response  of  the  linear  model  of  the  truss 
augmented  with  the  hnear  AJT  models  is  also  calculated  for  comparison.  The  results 
are  shown  in  figure  2.17.  The  results  show  a  very  close  agreement  in  the  overall 
shape  of  the  magnitude  plot,  especially  as  to  the  locations  of  the  poles  and  zeros. 
This  further  confirms  our  assumptions  made  when  adding  damping  to  the  model. 
The  phase  of  the  linear  model  drops  more  quickly  than  the  experimental  data-this  is 
due  to  the  model  of  the  AJTs  and  is  not  a  problem  with  the  linear  truss  FEM.  The 
bandwidth  of  the  AJTs  may  be  clearly  seen,  as  the  phase  and  magnitude  of  the  truss 
plummet  at  120  rad/sec  in  the  experimental  data,  a  phenomenon  reflected  by  the 
linear  model  which  rolls  off  at  about  40  dB/decade  and  seeks  360  degrees  of  phase 
lag. 

With  the  truss  and  actuator  models  thus  rigorously  verified,  we  may  now  turn  our 
attention  to  the  control  problem  of  vibration  suppression  and  expect  good  agreement 
between  theoretical  design  and  actual  experimental  results. 


recent  data.  That  explains  the  discrepancy  between  the  simulation  and  the  experimental  results. 


38 


frequency  (rad/sec) 


frequency  (rad/sec) 

Figure  2.17:  Open  loop  frequency  response 


Chapter  3 


CLASSICAL  CONTROL  TECHNIQUES  FOR 
VIBRATION  SUPPRESSION 


3.1  Overview 

Design  of  a  classical  controller  using  the  available  means  of  actuation  to  control  struc¬ 
tural  vibration  was  the  primary  goal  of  this  thesis.  Velocity  feedback  controllers  of 
the  type  used  at  the  Air  Force  Academy  are  presented  here,  complete  with  discussions 
of  the  theoretical  basis  of  using  velocity  control. 

3.2  Velocity  Feedback  to  the  AJTs 

Velocity  feedback  for  vibration  suppression  has  long  been  recognized  as  a  primary 
means  of  achieving  higher  damping  ratios  in  lightly  damped  systems.  Zimmerman 
applied  velocity  feedback  to  damp  the  transverse  (bending)  modes  of  a  0.654m  can¬ 
tilever  beam  with  a  reaction  mass  actuator  [23].  The  USAFA  truss  was  controlled 
using  velocity  feedback  to  both  RMAs  and  AJTs  [8].  It  is  worthwhile  to  explore  some 
of  the  guarantees  and  shortcomings  of  velocity  feedback,  and  to  report  test  results 
performed  on  the  UWAA  planar  truss. 

3.2.1  The  Theory  of  Velocity  Feedback 

Consider  a  spring-mass-damper  system  with  an  available  forcing  function  such  as  the 
reaction  mass  actuator  in  figure  2.5.  The  system  has  been  shown  to  be  modeled  by 
the  differential  equation: 


m„x  +  c^x  +  kaX  =  (3.1) 

Control  of  the  mass  is  easily  achieved  if  the  velocity  signal  is  present  for  feedback 
to  the  force  generator  (in  the  case  of  the  RMA,  feedback  to  the  coil’s  amplifier).  As 
shown  by  Inman  and  others  [10],  the  feedback  is  of  the  form  of  —KGx  where  G  is  the 


gain  constant  cf  the  RM  A  and  K  is  the  amplifier  gain.  Dropping  the  subscripts  from 
equation  3.1  and  letting  KG  =  g  the  differential  equation  for  the  system  undergoing 
feedback  is: 

mx  4-  ci  +  fcx  =  —gx  (3.2) 

Rewriting  equation  3.2  in  standard  second  order  form  [10]: 

X  4-  — — x  +  =  0  (3.3) 

m 

Equation  3.2  may  be  rewritten  as: 

X  4-  ^-i^x  4-  —  =  0  (3.4) 

m  m 

where  the  natural  frequency  is  a;„  =  yjkjm  and  the  damping  ratio  is  (^  =  (c  4- 
g)l2\lklm.  We  see  that  for  any  amplifier  gain  KG  ■=  g  negative  velocity  feedback 
will  result  in  a  stable  system  since  k,  c,  and  m  are  all  positive  values. 

Although  this  is  a  simple  example,  it  illustrates  the  advantage  of  velocity  feedback. 
For  a  reaction-less  actuator  without  any  lag  in  its  application  of  force,  any  velocity 
feedback  gain  will  serve  to  increase  the  stability  of  the  system  to  which  it  is  applied. 
This  is  the  main  impetus  behind  using  velocity  feedback  to  the  AJTs  for  vibration 
suppression. 

3.2.2  Approximate  Integrator  for  Velocity  Feedback 

To  achieve  velocity  feedback,  a  velocity  signal  is  obviously  required.  On  the  UW 
planar  truss  in  the  current  configuration  the  accelerometers  are  colocated  with  the 
AJTs  at  the  truss  tip  and  midpoints.  These  sensors  measure  the  lateral  acceleration 
of  the  truss  at  their  respective  locations.  To  achieve  a  velocity  signal,  the  first  guess 
would  be  to  try  to  integrate  the  acceleration  signal  directly.  There  are  problems  with 
this  set-up  as  any  slight  bias  in  the  system  while  the  planar  truss  is  in  the  unexcited, 
resting  state  would  be  integrated  until  a  sufficient  level  of  voltage  is  built  up,  tripping 
the  switching  transistor  circuit  of  the  AJTs  and  forcing  an  unwarranted,  destabilizing 
firing  of  the  AJT  associated  with  the  runaway  bias  signal.  Therefore,  some  feedback 
is  desirable  to  keep  unwanted  bias  from  triggering  the  system  at  rest. 
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The  USAFA  research  team  used  an  approximate  integrator  of  the  type  used  by 
the  Martin  Marietta  Astronautics  Group  [8].  The  approximate  integrator  is  modeled 
by  the  differential  equation: 


+  (leyi  +  Ci^Cyi  — 


(3.5) 


where  Cai  is  the  accelerometer  signal,  e„«  is  the  velocity  signal,  Q  is  a  user  constant,  and 
RC  is  a  user  specified  time  constant.  Taking  the  Laplace  transform  of  equation  3.5 


freq  (rad/sec) 

Figure  3.1:  Approximate  integrator  theoretical  frequency  response 


yields 


(3,6) 
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which  may  be  put  into  the  transfer  function  form: 


G(s) 


_ s/^ _ 


(3.7) 


The  frequency  response  of  the  transfer  function  G(s)  is  shown  in  figure  3.1  for 
constants  of  BC  =  1  and  fl  =  0.628.  The  resulting  frequency  response  is  very  close 
to  that  of  an  ideal  integrator’s  frequency  response  for  frequencies  above  10  rad/sec. 


Figure  3.2:  Approximate  integrator  analog  computer  implementation 

To  implement  equation  3.5  a  state  space  representation  is  more  convenient  to  use. 
The  states  may  be  defined  by  the  following  equations  [8]: 


evi 

where  e^ 
or  62 


-Qevi  +  fi  + 

— fiev, 

— fi  J  evidt  -f-  i.c. 


The  state  space  equations  which  result  are: 


eVi 

^2 


~n  n' 

evi 

-Q  0 

£2 

(3.8) 
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y 


+  0  ea, 


(3.9) 


Figure  3.2  in  an  analog  computer  diagram  of  equation  3.9  which  is  suitable  for  pro¬ 
gramming  the  GP-6  and  EAI-20.  The  values  for  f)  again  equals  0.628  and  the  constant 
RC  again  has  the  value  of  one.  It  is  important  to  note  that  the  integrators  and  sum¬ 
mers  of  an  analog  computer  invert  the  incoming  signal.  Figure  3.2  shows  a  negative 
(inverted)  velocity  signal  as  the  system  output  rather  than  the  positive  velocity  sig¬ 
nal  derived  in  the  state  space  equations.  This  is  because  it  is  convenient  to  allow 
the  velocity  signal  to  remain  negative  for  feedback  purposes.  If  a  positive  signal  was 
desired,  an  inverter  would  need  to  be  added  at  the  output. 


3.2.3  AJT  Control:  Non-Linear  Simulation  and  Experimental  Results 

The  classical  controller  for  the  AJT  control  circuit  is  simply  velocity  feedback  to  the 
AJT  control  circuitry.  As  mentioned  in  section  3.2.1,  negative  velocity  feedback  will 
always  stabilize  a  system  when  there  is  no  appreciable  lag  in  the  feedback  loop.  Un¬ 
fortunately,  there  is  significant  phase  lag  in  the  AJTs.  However,  for  the  low  frequency 
bending  modes,  the  AJTs  still  provide  a  significant  increase  in  damping. 


Table  3.1:  Non-linear  simulation  closed  loop  modal  data:  AJT  velocity  feedback 


frequency  (rad/sec) 

damping  ratio 

mode  1 

9.30 

0.06 

0.945 

mode  2 

55.80 

0.0345 

0.915 

mode  3 

146.52 

0.009 

0.000 

mode  4 

262.32 

0.009 

-0.250 

Using  non-linear  simulation  with  the  theoretical  non-linear  AJT  model  acting  upon 
the  truss  tip  and  midpoint  lateral  DOF  the  time  response  in  figure  3.3  were  obtained. 
The  truss  was  driven  in  same  manner  as  the  experimental  data,  a  sinusoidal  wave 
was  input  to  the  A  JT  circuit  at  the  different  modal  frequencies  of  the  linear  model. 
However,  instead  of  allowing  the  truss  to  reach  steady  state  excitation  as  is  the  Cctse 
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seconds 


Figure  3.3:  Non-linear  simulation  truss 


tip  responses.  AJT  velocity  feedback 


with  the  experimental  data,  the  sine  signal  was  switched  to  the  control  signal  when  the 
modal  excitation  of  the  linear  model  was  approximately  the  same  value  as  obtained 
in  the  experimental  testing.  The  modal  decay  is  seen  to  occur  at  a  much  faster  rate 
for  modes  one  and  two  than  in  the  open  loop  case  presented  in  figure  2.15.  This  is 
expected  due  to  the  AJT’s  lack  of  phase  lag  at  the  lower  modes.  A  small  oscillation 
is  still  apparent  in  the  second  mode  plot  as  the  AJTs  stop  firing  when  the  control 
signal  drops  into  the  deadband  voltage  cut-off  level.  Modes  three  and  four  are  not 
visibly  affected  by  the  addition  of  the  AJT  controllers. 

To  help  to  quantify  the  effects  of  the  controller,  damping  ratios  were  calculated 
for  each  mode  utilizing  the  logarithmic  decrement  method  of  section  2.5.2  and  are 
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Figure  3.4;  Non-linear/hardware  block  diagram  for  AJT  velocity  feedback  control 

displayed  in  table  3.1.  It  is  realized  that  no  linear  method  imposed  on  a  non-linear 
system  can  fully  quantify  the  damping  as  all  damping  is  assumed  to  be  viscous,  but 
it  is  still  valuable  for  comparison  to  the  open  loop  system.  The  log^o  of  the  ratio  of 
the  closed  loop  to  the  open  loop  value  is  displayed  in  table  3.1  in  favor  of  percentages 
because  of  the  large  increases  in  damping  achieved  (nearly  a  factor  of  10  for  both 
the  first  and  second  modes  in  non-linear  simulation).  The  damping  of  the  third  and 
fourth  modes  is  not  positively  affected,  in  fact  the  damping  decreases  slightly  in  the 
fourth  mode  compared  to  the  open  loop  system. 

The  experimental  truss  tip  responses  are  presented  in  figure  3.5,  and  show  a  slightly 
different  story  than  predicted  by  the  non-linear  simulation  of  AJT  velocity  feedback 
performance.  The  third  and  fourth  mode  have  AJT  misfires  which  are  destabilizing, 
indicating  that  the  AJT  velocity  feedback  controller  is  detrimental  at  higher  frequen¬ 
cies,  not  of  little  effect  as  predicted  by  the  non-linear  simulation.  The  first  and  second 


seconds  seconds 

Figure  3.5:  Experimental  truss  tip  responses:  AJT  velocity  feedback 

modes  behave  similarly  to  the  non-linear  simulation  predicted  results,  including  the 
small  oscillation  of  the  second  mode  when  the  control  signal  is  within  the  AJT  dead¬ 
band.  The  second  mode  oscillation  was  noticed  in  tj'\e  USAFA  truss  experiments,  but 
specific  mention  was  made  that  the  AJTs  never  destabilized  a  mode  [8].  It  seems 
that  either  slightly  greater  lags  in  our  AJTs  or  slightly  more  noise  in  our  control 
circuitry  must  have  been  present.  As  seen  in  table  3.2,  significant  damping  increases 
were  reported  for  first  and  second  mode  damping  ratios,  although  slightly  below  the 
performance  achieved  in  non-linear  simulation.  The  first  and  second  mode  damping 
increases  were  in  the  same  relative  proportion  at  about  7/lOths  of  a  factor  of  ten, 
just  as  in  non-linear  simulation  the  increases  were  in  the  same  relative  proportion  at 
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9/lOths  of  a  factor  of  10 

Table  3.2:  Experimental  closed  loop  modal  data:  truss  velocity  feedback  to  AJTs 


frequency  (rad/sec) 

damping  ratio 

mode  1 

9.30 

0.035 

0.712 

mode  2 

55.80 

0.025 

0.775 

mode  3 

146.52 

AJT  misfires 

na 

mode  4 

262.32 

AJT  misfires 

na 

Overall,  it  was  demonstrated  that  the  AJTs  were  effective  in  controlling  low  fre¬ 
quency  first  and  second  mode  disturbances.  However,  the  failure  to  effectively  dampen 
the  vibrations  of  the  third  and  fourth  modes,  and  possibly  excite  them,  indicates  that 
the  use  of  velocity  feedback  to  the  AJTs  alone  is  an  inadequate  solution  to  the  vibra¬ 
tion  suppression  problem. 

3.3  RMAs  for  Vibration  Suppression 

The  reaction  mass  actuators  are  linear,  space-realizable  forms  of  actuation  ideally 
suited  to  control  lateral  vibration  modes  of  the  planar  truss.  Unlike  the  AJT  thrusters, 
only  electrical  energy  is  necessary  to  power  these  actuators,  not  the  expulsion  of  a 
finite  supply  of  gas  as  is  the  case  with  the  AJTs.  However,  the  amount  of  force  the 
RMAs  may  exert  is  somewhat  limited,  and  they  are  limited  in  effectiveness  by  the 
their  non-tuneable  natural  frequencies  as  we  shall  see.  The  RMAs  are  located  at  the 
truss  midpoint  and  tip,  locations  which  might  be  adjusted  later  to  suppress  different 
modes.  As  the  RMAs  are  set  up  in  the  current  configuration,  the  midpoint  RMA  is 
at  an  anti-node  of  the  third  lateral  truss  bending  mode,  which  limits  its  effectiveness 
to  control  that  particular  mode. 

A  discussion  of  the  different  levels  of  classical  control  is  given  in  this  section. 
The  use  of  the  RMAs  in  the  open  loop  system  as  passive  dampers  is  described  in 
section  3.3.1.  Feedback  of  only  the  truss’  velocity  signals  as  performed  by  Zimmer¬ 
man  [23]  is  discussed  in  section  3.3.2.  The  use  of  relative  velocities  between  the 
truss’  inertial  velocities  and  the  RMAs’  proof  mass  inertial  velocities  is  presented  in 
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section  3.3.3  and  builds  upon  the  separate  works  of  Hallauer  [8],  Miller  and  Craw¬ 
ley  [15],  and  Inman  [10].  Finally,  an  approach  proposed  by  Inman  which  utilizes 
separate  gains  on  the  proof  mass  velocities  and  the  truss  velocities  is  applied  to  yield 
a  very  attractive  classical  velocity  control  feedback  system.  These  results  are  reported 
in  section  3.3.4. 


3.3.1  RMA  as  a  Passive  Damper 

The  use  of  reaction  mass  actuators  as  vibration  absorbers  is  not  a  new  concept. 
The  simplicity  of  merely  adding  a  spring-mass-damper  to  a  dynamic  system  to  re¬ 
duce  vibrations  is  very  attractive,  and  the  problem  has  been  thoroughly  explored. 
Applications  to  large  space  structures,  such  as  the  optimal  design  of  a  vibration  ab¬ 
sorbers  for  multi-degree  of  freedom  (MDOF)  systems  such  as  beams  and  trusses, 
have  been  the  subject  of  recent  literature  [11,  15].  The  results  give  general  guidelines 
for  the  application  of  vibration  absorbers  to  systems  similar  to  the  UWAA  planar 
truss.  Although  Miller  and  Crawley  [15]  used  quadratic  optimization  techniques  to 
arrive  at  a  closed  form  solution  to  the  passive  damping  problem,  their  results  may 
be  summarized  as  follows.  A  vibration  absorber  may  be  added  to  a  MDOF  system 
to  reduce  the  response  of  structural  modes.  Tuning  the  vibration  absorber  to  control 
the  vibration  of  a  system  is  performed  by  selecting  the  lowest  frequency  mode  of  in¬ 
terest  and  matching  the  natural  frequency  of  the  absorber  to  that  natural  frequency. 
The  vibration  absorber  is  placed  at  point  of  maximum  modal  excitation  for  best  re¬ 
sults.  The  frequency  of  the  vibration  absorber  will  have  to  be  slightly  lessened  as  its 
mass  approaches  the  equivalent  modal  mass  for  the  structural  mode  of  interest.  The 
damping  ratio  of  the  vibration  absorber  is  then  increased  until  the  desired  structural 
modal  damping  is  achieved.  Many  vibration  absorbers  may  be  individually  tuned 
and  placed  strategically  along  the  structure  to  attenuate  multiple  modal  responses, 
realizing  that  slight  residual  damping  will  occur  in  the  modes  above  the  tuned  fre¬ 
quency,  while  modes  below  the  tuned  frequency  will  be  slightly  adversely  affected. 

For  the  UWAA  planar  truss  with  the  RMAs  acting  at  the  lateral  truss  tip  and 
endpoint  degrees  of  freedom,  the  results  of  Miller  and  Crawley  may  be  used  to  predict 
the  effect  of  the  passive  damping  of  the  RMAs  on  the  planar  truss.  The  RMAs  are 
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Table  3.3:  Open  loop  linear  model  modal  data:  RMAs  as  passive  dampers 


frequency  (rad/sec) 

damping  ratio 

mode  1 

9.19 

0.0067 

-.006 

mode  2 

57.95 

0.0190 

0.656 

mode  3 

146.58 

0.0098 

0.037 

mode  4 

260.30 

0.0164 

0.011 

RMA 

40.64 

0.0635 

-0.073 

RMA 

43.19 

0.0776 

0.000 

acting  on  positions  of  maximum  displacement  for  the  second  mode  (see  figure  2.13). 
Furthermore,  the  natural  frequencies  of  the  RMAs  occur  at  41.12  rad/sec  for  the  tip 
RMA  and  42.96  rad/sec  for  the  midpoint  RMA,  both  of  which  are  about  10  rad/sec 
below  the  second  mode  but  well  above  the  first  bending  mode.  Therefore,  we  expect 
the  RMAs  will  add  significant  damping  to  the  second  bending  mode  of  the  truss  but 
below  their  true  attainable  capability  if  they  were  to  have  their  natural  frequencies 
lowered  by  the  addition  of  some  more  mass  (since  the  spring  constants  of  the  RMAs 
are  fixed).  Since  the  hardware  is  already  assembled,  the  masses  will  not  be  changed 
for  this  experiment,  but  leaves  room  for  a  study  to  be  performed  where  the  mass  of 
the  RMAs  is  varied. 

The  linear  model  for  the  RMA  system  is  formed  using  the  methods  of  sections  2.3.1 
and  2.4.2.  It  is  a  twelfth  order  system  of  the  form 

X  =  Aux  +  B12U  (3.10) 


The  numerical  values  which  constitute  Au  and  Bu  are  listed  in  section  A.2. 

As  we  see  in  table  3.3,  the  predictions  that  the  damping  ratio  of  the  second  mode 
would  significantly  increcise  are  correct  as  a  0.656  factor  of  ten  increase  in  second  mode 
damping  occurred  in  the  linear  model.  The  goal  of  10%  damping  in  each  structural 
mode  is  not  achieved  with  the  addition  of  passive  damping  alone,  although  damping 
of  the  second  mode  is  increased  substantially.  The  next  step  is  to  add  active  control 
to  supplement  the  passive  damping. 
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3.3.2  Straight  Velocity  Feedback 

This  section  begins  the  discussions  on  classical  control  of  the  truss  through  the  use 
of  the  RMAs.  All  feedback  gains  in  this  and  following  sections  assume  the  use  of  the 
a  scaled  accelerometer  output  equal  to  38.8  volts  per  g,  or  approximately  0.1  volts 
per  1  ■^.  This  is  easily  achieved  in  the  hardware  as  the  amplifier  gain  on  the  RMA 
amplifier  is  marked  in  units  ranging  from  0  to  10.  Additionally,  negative  feedback  is 
assumed  unless  a  negative  gain  is  specified,  which  indicates  positive  feedback. 


truss  velocity  feedback  relative  velocity  feedback 


Real  Real 

Figure  3.6;  Comparison  of  root  loci:  kt  =  km  =  0.0:0.1:10 


The  use  of  velocity  feedback  from  the  truss  stations  without  regard  to  the  inertial 
or  relative  velocities  of  the  RMA  proof  mass  is  similar  to  the  work  of  Zimmerman  [23] 
where  control  of  a  cantilever  beam  was  performed  using  an  RMA  tuned  to  the  first 
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bending  mode  of  the  beam.  Simple  proportional  feedback  of  the  beam  tip  velocity  to 
the  RMA  amplifier  yields  a  feedback  law  of  the  form: 

fb(tip)  —  ^tQt 

The  beam  was  then  only  excited  at  the  first  mode  at  amplitudes  well  within  the 
stroke  limits  of  the  RMA.  Two  problems  were  apparent.  First,  there  was  no  way  to 
predict  how  different  feedback  gains  would  affect  the  throw  of  the  proof  mass  since 
actuator  dynamics  were  ignored  in  the  modeling  as  well  as  in  the  implementation. 
Second,  small  levels  of  vibration  persisted  well  after  the  most  of  the  modal  energy 
was  dissipated.  Although  it  was  felt  that  the  D/A  converters  were  responsible  for  the 
problem,  the  fact  that  straight  velocity  feedback  is  risky  and  assumes  a  well  behaved 
system  to  begin  with  may  have  well  compounded  the  D/A  problem. 

Table  3.4:  Feedback  gains  on  truss  velocities 


velocity 

feedback  gain 

value 

4t 

kf 

2.0 

9m 

k 

2.0 

Figure  3.6  illustrates  the  disadvantage  of  the  using  straight  truss  station  velocity 
feedback  when  compared  to  feedback  which  takes  into  account  the  dynamics  of  the 
proof  mass  on  board  the  RMA.  The  feedback  control  law  was  of  the  form 

fb(tip)  — 

fb(mid)  —  (3.12) 

The  root  locus  for  feedback  gains  of  0  to  10  on  the  truss  tip  and  midpoint  velocities 
only  indicates  that  the  implementation  of  the  feedback  law  results  in  significant  stabi¬ 
lization  of  only  the  second  mode.  The  gains  in  table  3.4  are  chosen  where  the  second 
mode  damping  is  the  highest.  Although  the  damping  of  mode  1  is  also  improved,  the 
trade  off  versus  relative  velocity  feedback  is  the  complete  lack  of  authority  over  the 
higher  modes  of  the  system. 
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Table  3.5:  Closed  loop  linear  model  modal  data:  truss  velocity  feedback 


frequency  (rad/sec) 

damping  ratio 

mode  1 

9.19 

0.0100 

0.167 

mode  2 

56.85 

0.0298 

0.850 

mode  3 

146.57 

0.0098 

0.037 

mode  4 

260.28 

0.0167 

0.019 

RMA 

41.43 

0.2972 

0.598 

RMA 

43.19 

0.31466 

0.608 

The  use  of  straight  velocity  feedback  does  have  one  very  significant  advantage:  ease 
of  application.  The  number  of  sensors  required  is  one-half  the  number  required  for 
relative  velocity  feedback.  However,  given  that  the  UWAA  planar  truss  is  representa¬ 
tive  of  a  high  performance  requirement  large  space  structure,  that  advantage  does  not 
outweigh  the  disadvantage  of  overall  poor  performance.  It  was  decided  not  to  imple¬ 
ment  the  velocity  controller  in  the  hardware  because  of  the  difficulty  in  determining 
a  performance  increase  over  the  passively  damped  system  of  section  3.3.1. 

3.3.3  Relative  Velocity  Feedback 

The  use  of  relative  velocity  feedback  was  used  effectively  by  Hallauer  to  suppress 
bending  mode  vibrations  [8|.  The  feedback  control  law  is  of  the  form 

fb(tip)  —  ktQta 

fb{mid)  —  km(Jrna  (3.13) 

where  and  q,na  are  the  relative  velocities  between  the  tip  actuator  mass  and  tip 
truss  station  and  the  mid  acutuator  mass  and  mid  truss  station,  respectively. 

As  seen  in  figure  3.6  the  feedback  of  relative  velocities  to  the  RMAs  has  a  much 
greater  effect  on  the  closed  loop  eigenvalues  of  the  truss  than  does  the  straight  velocity 
feedback  of  the  previous  section.  The  feedback  gains  for  the  truss  and  midpoint 
RMAs  were  adjusted  in  equal  increments  to  provide  the  levels  of  damping  which 
are  shown  graphically  in  figure  3.7.  The  first,  third,  and  fourth  bending  modes 


damping  ratio  damping  ratio 

Figure  3.7;  Relative  velocity  feedback  to  the  RMAs:  kt  =  km  =  0:0.4:10 

all  respond  favorably  to  increases  in  feedback  gain  within  the  range  of  the  RMA 
amplifiers.  Increases  in  gain  beyond  the  range  of  the  amplifiers  would  easily  be 
obtained  by  adding  additional  gains  in  the  feedback  loops  using  the  analog  computers, 
but  high  gain  systems  are  generally  undesirable  due  to  the  penalty  of  increased  noise 
in  the  feedback  loop  and  were  not  considered  here.  For  the  second  mode,  however, 
increasing  the  gain  adversely  affected  stability  after  a  value  of  0.2  was  reached.  Similar 
behavior  occurs  for  the  third  and  fourth  modes  at  higher  levels  of  gain  avoided  here. 
This  led  Hallauer  to  correctly  state  that  the  best  performance  using  relative  velocity 
feedback  would  be  obtained  if  the  gains  were  implemented  using  band-pass  filters 
which  would  allow  for  the  best  possible  gain  individually  selected  for  each  mode  [8]. 
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frequency  (rad/sec) 


frequency  (rad/sec) 

Figure  3.8:  Linear  model  frequency  response  for  mode-2  suppression:  kt  =  km  =  2.0 


This  approach,  however,  is  very  difficult  to  implement  in  the  hardware. 

The  ability  of  the  RMAs  to  suppress  vibration  when  gains  are  chosen  for  a  specific 
mode  is  shown  in  the  frequency  responses  of  figure  3.8  and  figure  3.9.  The  data  in 
figure  3.8  was  obtained  by  using  the  linear  model  of  the  AJTs  to  input  force  to  the 
RMA-augmented  truss  model  with  the  feedback  gains  on  the  relative  velocities  set  to 
2.0  in  the  closed  loop  case.  Comparison  of  the  open  and  closed  loop  responses  shows 
6.3  dB  of  vibration  suppression  occurred  at  the  second  mode.  A  similar  approach  was 
applied  in  the  case  of  the  experimental  hardware  results  shown  in  figure  3.9  where 
suppression  of  the  second  mode  due  to  AJT  tip  excitation  is  8.58  dB.  The  closed  loop 
gains  for  the  experimental  case  were  found  by  exciting  the  truss  at  the  second  mode 
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frequency  (rad/sec) 


frequency  with  the  AJTs  and  then  adjusting  the  amplifier  gains  and  measuring  the 
real-time  acceleration  output  of  the  truss  tip  and  midpoint  accelerations  until  the 
accelerations  were  at  a  minimum.  The  experimental  data  confirms  the  linear  model’s 
utility  in  predicting  second  mode  vibration  suppression. 

As  previously  stated,  to  design  a  band-pass  filter  for  individual  modes  makes  a 
difficult  task  for  hardware  implementation.  Instead,  it  is  desirable  to  have  one  set  of 
gains  which  may  be  implemented  which  might  realize  most  of  the  performance  of  the 
band-pass  filter  designs.  Taking  the  approach  that  the  midpoint-mounted  RMA  has 
little  effect  on  the  third  bending  mode  of  the  truss  due  to  its  location  at  a  third  mode 
node,  the  relative  velocity  feedback  gain  was  set  to  a  value  to  dampen  the  second 
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Table  3.6:  Feedback  gains  on  relative  velocities 


velocity 

feedback  gain 

value 

Qa  ~  Qt 

kt 

6.0 

3a  9m 

1.9 

mode  of  the  truss.  Then,  to  increase  the  damping  of  the  third  and  fourth  modes,  the 
tip- mounted  RMA  had  its  gain  increased  to  values  more  in  line  with  damping  the 
higher  modes.  The  final  gains  chosen  are  shown  in  table  3.6.  The  modal  damping  is 
shown  in  table  3.7  for  the  linear  model  results. 

Table  3.7:  Closed  loop  linear  model  modal  data:  relative  velocity  feedback 


frequency  (rad/ sec) 

damping  ratio 

mode  1 

9.20 

0.013 

mode  2 

55.12 

0.943 

mode  3 

146.04 

0.0156 

0.239 

mode  4 

260.00 

0.0205 

0.108 

RMA 

43.04 

0.2904 

0.586 

RMA 

43.08 

0.7764 

1.000 

This  is  a  suitable  compromise  as  the  damping  iS'  significantly  increased  for  modes 
2,  3,  and  4.  The  damping  of  the  second  mode  is  still  nearly  a  factor  of  ten  greater  than 
the  open  loop  case,  while  the  damping  of  the  third  and  fourth  modes  are  increased 
by  a  smaller  factor.  This  makes  sense  as  the  third  and  fourth  modal  frequencies  are 
’further’  away  from  the  ideal  controlling  frequency  of  the  RMAs  which  is  just  below 
the  natural  frequency  of  the  second  mode. 

3.3.4  Separate  Velocity  Feedback  using  RMA  and  Truss  Inertial  Velocities 

Inman  stated  in  the  conclusions  of  his  paper  on  control/structure  interaction  [10],  that 
a  better  approach  for  using  the  reaction  mass  actuator  would  be  to  use  a  separate 
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open  loop  (x),  closed  loop  rel.  vel.  only  (+),  k  =  0.0:0.05:0.7 
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Figure  3.10:  Root  locus  for  separate  feedback  g,a.ins  on  inertial  velocities 

feedback  path  for  the  inertial  velocities  of  the  proof  mass  and  the  corresponding  truss 
station.  The  feedback  law  for  the  truss  tip  would  take  the  form  of: 

fb  =  -kttqt  -  ktqa  (3.14) 

where  qt  is  the  truss  tip  inertial  velocity  and  (ja  and  proof  mass  inertial  velocity.  A  root 
locus  of  showing  the  effects  of  separating  the  feedback  paths  is  shown  in  figure  3.10. 
The  root  locus  is  obtained  by  beginning  with  the  closed  loop  system  using  relative 
velocity  feedback  as  presented  in  section  3.3.3  (marked  as  Then,  the  inertial 

truss  velocity  is  added  to  the  relative  velocity  as  the  root  locus  gain  is  increased: 
(ja  —  k(jf  The  first  mode  is  destabilized  by  this  type  of  feedback,  so  a  limit  of  0.7  was 


6-5-4-3-2-10  1 

Real 
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chosen  for  k  to  preserve  the  stability  of  the  first  mode.  Since  the  ultimate  goal  of 
this  research  was  to  improve  damping  using  both  the  RMAs  and  AJTs,  the  increase 
in  the  higher  mode  damping  where  the  AJTs  are  ineffective  is  seen  as  a  worthwhile 
trade-off  for  the  decrease  in  the  first  mode  damping  where  the  AJT  performance  is 
most  effective. 


seconds  seconds 

Figure  3. IT.  Time  responses  for  second  mode;  separate  inertial  velocity  feedback 

The  feedback  gains  on  the  inertial  velocities  chosen  for  the  final  design  are  pre¬ 
sented  in  table  3.8.  The  resulting  modal  damping  from  the  linear  model  is  shown 
in  table  3.9.  The  damping  of  the  second,  third,  and  fourth  modes  is  approximately 
50  %  better  than  the  damping  achieved  through  the  use  of  relative  velocity  feedback 
alone  (compare  to  table  3.7).  A  time  response  for  the  linear  model  second  mode 
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Table  3.8:  Feedback  gains  on  proof  mass  and  truss  velocities 


velocity 

feedback  gain 

value 

qa{tip) 

kt 

6.0 

qt 

kit 

-1.2 

qaimid) 

1.9 

qm 

-0.57 

Table  3.9:  Closed  loop  linear  model  modal  data:  separate  inertial  velocity  feedback 


frequency  (rad/ sec) 

damping  ratio 

mode  1 

9.18 

0.00146 

-0.668 

mode  2 

54.24 

0.0446 

1.026 

mode  3 

145.72 

0.0192 

0.329 

mode  4 

259.82 

0.0230 

0.158 

RMA 

43.29 

0.2882 

0.585 

RMA 

43.76 

0.8856 

1.057 

vibration  suppression  is  shown  in  figure  3.11.  The  figure  shows  the  open  loop  free 
decay  and  the  closed  loop  time  response  for  the  truss  tip  acceleration.  It  also  shows 
the  reaction  mass  relative  displacement  for  the  tip  and  midpoint  RMAs.  These  plots 
show  that  the  reaction  masses  are  both  well  within  their  limited  throw  excursions  of 
±  0.1  in.  In  fact,  it  would  take  second  mode  vibrations  well  in  excess  of  100  to 
cause  the  reaction  masses  to  travel  beyond  their  throw  limitations.  The  second  mode 
plot  was  included  because  the  amplitude  of  travel  of  the  reaction  mass  was  greater  at 
this  mode  than  the  higher  frequency  third  and  fourth  modes.  The  first  mode  is  not 
examined  because  the  RMAs  do  not  aid  the  damping  of  that  mode,  although  the  tip 
RMA  will  reach  the  limits  of  its  throw  excursion  for  first  mode  vibrations  in  excess 
of  50  measured  at  the  tip  accelerometer. 
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3.4  Hybrid  Control:  AJTs  and  RMAs  Working  Together  to  Suppress  Vibration 

Hybrid  actuation  was  the  term  used  by  Hallauer  [8]  to  describe  the  use  of  both  the 
AJTs  and  the  RMAs  to  suppress  vibration  of  the  planar  truss  bending  modes.  To 
employ  the  two  different  actuators  in  the  non-linear  simulation  and  experimental 
hardware,  it  is  first  necessary  to  note  that  both  of  the  controllers  were  designed  to 
work  independently  of  each  other.  No  effort  was  made  to  make  the  controllers  work 
together  except  the  reasoning  in  that  the  AJTs  would  dampen  the  low  frequency 
vibrations  while  the  RMAs  would  suppress  the  high  frequency  vibrations.  This  control 
effectively  ’piggy-backs’  the  non-linear  AJT  controllers  on  top  of  the  RMA-augmented 
linear  closed  loop  system.  Effectiveness  of  this  system  will  allow  the  design  of  an 
optimal  controller  for  the  linear  RMA-augmented  truss  without  regard  to  the  AJTs, 
which  then  may  be  piggy-backed  on  board  the  optimally  controlled  truss. 

Table  3.10:  Non-linear  simulation  closed  loop  modal  data:  hybrid  actuation 


frequency  (rad/sec) 

damping  ratio 

UovenlooJ 

mode  1 

9.18 

0.0904 

1.12 

mode  2 

54.24 

0.0914 

1.33 

mode  3 

145.72 

0.0170 

0.276 

mode  4 

259.82 

0.0212 

0.122 

The  non-linear  simulation  conducted  to  test  the  effectiveness  of  hybrid  actuation 
was  performed  by  using  the  AJT  velocity  controller  discussed  in  section  3.2.3  and 
the  RMA  controller  from  section  3.3.4.  Figure  3.13  shows  the  set  up  put  into  the 
non-linear  simulation  program.  The  RMA-controlled  truss  was  driven  with  the  AJTs 
at  the  natural  frequencies  until  a  level  of  excitation  which  was  representative  of  the 
experimental  excitations  weis  reached,  and  then  the  AJTs  control  circuit  was  switched 
on.  This  allowed  both  the  AJTs  and  the  RMAs  to  work  together  to  control  the  truss 
model.  The  non-linear  simulation  of  the  truss  tip  acceleration  time  responses  for  the 
four  bending  modes  is  presented  in  figure  3.12. 

From  the  data  presented  in  figure  3.12,  logarithmic  decrements  were  measured  and 
damping  ratios  calculated.  These  values  are  presented  in  table  3.10.  The  damping  of 


61 


0  10  20  0  2  4 


seconds  seconds 


seconds  seconds 

Figure  3.12;  Non-linear  simulation  of  hybrid  controhtruss  tip  responses 

the  modes  increases  as  is  expected.  In  general,  the  second,  third,  and  fourth  modes 
had  levels  of  damping  comparable  to  the  RMA-only  controlled  system  in  section  3.3.4, 
while  the  first  mode  damping  was  clearly  more  influenced  by  the  AJTs.  The  level  of 
damping  appears  to  increase  significantly  over  the  AJT-only  control  of  section  3.2.3, 
but  this  may  in  part  be  due  to  the  lower  amplitudes  of  the  hybrid  testing.  The 
damping  in  the  first  two  modes  nearly  meets  the  goal  of  ten  percent  damping  in 
future  large  space  structures.  The  next  step  is  to  perform  the  experimental  tests. 
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Figure  3.13:  Non-linear  block  diagram  for  hybrid  actuation  simulation 
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3.5  Experimental  Results 

3.5.1  Experimental  Open  Loop  Testing  of  the  RMA-Augmented  Truss 
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Figure  3.14:  Experimental  open  loop  RMA-augmented  truss  tip  responses 


The  experimental  testing  was  conducted  close  to  the  same  manner  as  the  non-linear 
simulation,  but  the  RMA-augmented  truss  was  excited  with  the  AJTs  without  any  of 
the  RMA  control  circuitry  on.  Open  loop  RMA-augmented  truss  plots  are  included 
from  the  day  of  the  experimental  testing  to  give  a  basis  of  comparison  between  the 
passive  damping  effects  of  the  RMAs  and  the  closed  loop  time  responses.  The  RMA- 
augmented  open  loop  plots  are  presented  in  figure  3.14.  The  displacement  differences 
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between  the  open  loop  time  responses  in  figures  2.15  and  2.16  are  due  to  two  reasons. 
One,  the  RMAs  prevent  the  full  level  of  truss  excitement  due  to  their  passive  damping 
effects.  Two,  the  testing  environmental  conditions  differed  between  the  two  test  days. 
It  was  cooler  in  temperature  when  the  RMA  tests  were  conducted  compared  to  the 
open  loop  tests  of  chapter  two. 

Table  3.11;  Experimental  open  loop  modal  data:  RMA-augmented  truss 


frequency  (rad/sec) 

damping  ratio 

mode  1 

9.18 

0.0120 

0.246 

mode  2 

54.24 

0.0156 

0.570 

mode  3 

145.72 

0.0085 

-0.025 

mode  4 

259.82 

0.0128 

-0.097 

The  damping  ratio  of  the  open  loop  RMA-augmented  truss  increase  at  the  second 
mode  of  the  truss  as  expected.  The  increase  in  damping  of  the  first  mode  is  not 
substantial  when  one  takes  into  consideration  the  change  in  damping  that  occured 
when  the  truss  was  retested  as  noted  in  section  2.5.2.  The  damping  in  mode  one 
is  identical  to  later  tests  on  that  mode  made  after  the  baseline  open  loop  damping 
was  established.  Therefore,  the  truss  behaves  as  expected  for  the  first  mode  with 
the  addition  of  the  RMA  passive  dampers.  The  change  in  damping  of  the  third  and 
fourth  mode  is  slight  as  predicted  by  the  modeling  performed  in  section  3.3.1. 

3.5.2  Experimental  Closed  Loop  Testing  of  the  RMA  Controllers 

The  next  step  in  the  experimental  testing  was  to  conduct  a  test  of  the  RMA  controller 
acting  independently.  A  problem  surfaced  when  amplifier  gains  were  set  at  levels 
above  7.0.  Although  this  problem  surfaced  when  the  feedback  loops  were  connected, 
open  loop  tests  confirmed  that  this  problem  was  isolated  to  the  amplifier  circuitry.  In 
figure  .3.15,  a  time  response  of  the  tip  actuator  operating  at  the  second  mode  of  the 
truss  is  displayed.  The  figure  shows  the  RMA  acceleration  compared  to  the  truss  tip 
acceleration.  For  the  first  second,  the  system  is  open-loop  with  the  gain  set  to  zero.  At 
about  one  second,  the  RMA’s  amplifier  is  switched  on  to  a  gain  of  seven  acting  on  the 
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relative  velocity  between  the  acutator  mass  and  the  truss  station.  Normal  attenuation 
is  recorded.  At  about  two  seconds  on  the  plot,  the  AJTs  providing  the  second  mode 
disturbance  to  the  truss  are  switched  off.  This  results  in  the  truss  station  acceleration 
decaying  to  zero,  but  an  actuator/instrumentation  instability  arises.  This  causes  the 
RMA  proof  mass  to  violently  swing  back  and  fourth,  even  when  the  steady  state  truss 
acceleration  is  zero.  A  problem  thought  to  be  similar  to  this  was  noticed  by  Miller 
and  Crawley  in  their  experiments  with  RMAs  [15]. 


time  (seconds) 


Figure  3.15:  RMA  actuator/instrumentation  instability 

Because  the  actuator/ instrumentation  instability  limited  the  gain  adjustments  pos¬ 
sible,  the  RMA  controller  for  experimental  implementation  was  chosen  to  be  the  ve¬ 
locity  feedback  controller  for  mode-2  suppression  presented  in  section  3.3.3.  The  tip 


seconds 


seconds 


Figure  3.16;  Experimental  RMA-relative  velocity  feedback  only-time  responses 

and  midpoint  amplifier  gains  were  set  to  values  of  3.5  for  best  experimental  results 
for  mode-2  damping. 

The  time  responses  of  the  experimental  RMA  controlled  truss  system  are  presented 
in  figure  3.16.  The  time  responses  show  the  RMAs  working  effectively  to  actively 
suppress  second  bending  mode  vibrations.  The  damping  ratios  were  calculated  and 
are  presented  in  table  3.12.  The  first  mode  was  nearly  unchanged  from  the  open  loop 
RMA-augmented  truss  system  in  section  3.5.1.  The  third  and  fourth  modes  both  have 
their  damping  increased  compared  to  the  open  loop  RMA-augmented  truss  as  was 
predicted  in  section  3.3.3.  The  second  mode  damping  is  most  noticeably  increased, 
as  was  indicated  by  the  modal  time  responses  in  figure  3.16. 
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Table  3.12:  Experimental  closed  loop  modal  data:  relative  velocity  feedback  to  RMAs 


frequency  (rad/sec) 

damping  ratio 

mode  1 

9.18 

0.0128 

0.275 

mode  2 

54.24 

0.0308 

0.865 

mode  3 

145.72 

0.0140 

0.192 

mode  4 

259.82 

0.0155 

-0.014 

3.5.3  Experimental  Hybrid  Control 

Implementation  of  the  experimental  hybrid  control  system  involved  using  the  AJT 
controller  as  reported  in  section  3.2.3  and  the  RMA  controller  as  reported  in  sec¬ 
tion  3.5.2.  A  schematic  of  the  controller  is  presented  in  figure  3.18.  The  dc  offset  to 
the  approximate  integrators  was  necessary  to  counter  slight  bias  in  the  accelerometer 
signals  resulting  from  the  measurement  of  a  small  component  of  gravity. 

Time  responses  were  obtained  by  driving  the  truss  at  the  respective  modal  frequen¬ 
cies  and  then  switching  both  the  RMA  control  circuts  and  the  AJT  control  circuits  on 
at  the  same  time.  The  modal  time  responses  are  shown  in  figure  3.17.  As  predicted 
by  the  non-linear  simulation  reported  in  section  3.4  and  by  the  research  performed  by 
Hallauer  on  the  USAFA  truss  [8],  the  damping  of  the  first,  second,  and  third  bending 
modes  is  increased.  The  AJTs  and  RMAs  work  particularly  well  at  damping  mode-2 
vibration.  The  AJT  still  misfires  during  mode-4  excitation,  but  the  RMAs  seemed  to 
have  kept  the  AJTs  from  misfiring  during  mode-3  vibration  suppression.  The  RMAs 
thus  appear  to  assist  the  AJTs  at  the  higher  modes.  The  fourth  mode  appears  to  be 
positively  affected  by  the  addition  of  the  RMAs  to  the  control  scheme  as  well  com¬ 
pared  to  the  AJTs  acting  alone.  Also,  the  mode-1  damping  appears  cleaner  when 
compared  to  the  AJT-only  control  in  figure  3.5.  The  RMAs  probably  damp  out  the 
high  frequency  noise  caused  by  the  bang-bang  control  of  the  AJTs.  The  damping 
ratios  are  shown  in  figure  3.13. 

The  damping  ratios  were  higher  than  the  levels  of  damping  achieved  by  either 
controller  acting  independently  except  for  the  fourth  mode  which  is  adversely  affected 
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seconds 


Figure  3.17:  Experimental  hybrid  control  truss  tip  time  responses 


by  misfirings  of  the  AJTs. 


3.6  Classical  Control  Summary 

The  use  of  classical  control  techniques  is  effective  for  vibration  suppression  of  the 
planar  bending  modes  of  the  UWAA  truss.  The  levels  of  damping  achieved  using  the 
techniques  presented  in  this  study  fall  short  of  achieving  the  goal  of  10%  of  modal 
damping  desired  in  large  space  structures,  but  still  increased  the  damping  of  mode-1 
by  nearly  a  factor  of  10,  the  mode-2  damping  by  more  than  a  factor  of  ten,  and 
quadruples  the  level  of  damping  of  the  third  mode. 

Furthermore,  the  results  of  this  chapter  show  the  validity  of  the  modeling  of  chapter 
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Table  3.13:  Experimental  closed  loop  modal  data:  hybrid  actuation 


frequency  (rad/ sec) 

damping  ratio 

Uovcnl<^^ 

mode  1 

9.18 

0.0555 

0.912 

mode  2 

54.24 

0.0900 

1.33 

mode  3 

145.72 

0.0327 

0.560 

mode  4 

259.82 

AJT  misfires 

na 

two.  The  non-linear  simulation  was  at  all  times  consistent  with  the  hardware  with 
two  exceptions:  the  failure  to  predict  the  aciuator/instn*mentation  instability  in  the 
RMAs  and  the  misfiring  of  the  AJTs  when  controlling  the  third  and  fourth  modes  of 
the  truss.  The  performance  of  the  simulation  was  in  general  slightl3'  better  than  the 
experimental  results.  This  may  be  attributed  to  non-linearities  not  accounted  for  in 
the  modeling  of  the  truss  and  the  actuators  and  ignoring  the  effects  of  noise  in  the 
system. 

It  appears  that  the  goal  of  10%  modal  damping  is  within  the  reach  of  given  actua¬ 
tors  using  classical  control  techniques  for  the  first,  second,  and  third  modes.  Tuning 
of  the  natural  frequencies  of  the  RMAs,  addition  of  filters  to  the  AJT  control  circuit, 
and  the  choice  of  acutator  location  could  make  the  10%  damping  goal  a  reality. 
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Figure  3.18:  Hardware  block  diagram  for  hybrid  actuation 


Chapter  4 

OPTIMAL  CONTROL  STRATEGIES 


4.1  Overview 

Chapter  three  reported  the  design  and  testing  of  a  classical  multiple-input,  multiple- 
output  hybrid  actuation  controller  which  delivered  appreciable  increases  in  damping 
of  the  lateral  bending  modes  of  the  UWAA  planar  truss.  In  this  chapter  the  applica¬ 
tion  of  linear  quadratic  Gaussian  techniques  to  design  a  full  state  feedback  controller 
for  the  reaction  mass  actuators  is  discussed,  along  with  the  non-linear  simulation  of 
a  state  feedback  RMA  controller  working  together  with  the  air  jet  thruster’s  velocity 
controller  previously  discussed  in  chapter  three. 

4.2  LQG  Design  Problem 

The  solution  of  the  linear  quadratic  Gaussian  problem  is  well  understood[13].  Given 
the  plant  model  in  state-space  form: 


X  =  Ax  -t-  Bu  -t-  Ftc 

,V  =  Cx -f  n  (4.1) 

where  w  and  v  are  zero-mean  Gaussian  stochastic  processes  (commonly  referred  to 
as  white  noises)  assumed  to  be  uncorrelated  such  that: 

EficiT}  =  0  (4.2) 

The  problem  is  tlieri  to  form  a  feedback  control  law  which  minimizes  the  cost  func¬ 
tion  [13]: 

./  =  -  lim  E{  f  -f-  Iiu)dt} 

2  ^  Jo 


(4.3) 
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and  2  is  a  combination  of  the  states  given  by: 

2  =  Mx  (4.4) 

Q  is  a  positive  semi-definite  weighting  (penalty)  matrix  on  the  linear  combinations  of 
the  state  variables  and  7?  is  a  positive  definite  weighting  matrix  on  the  control  inputs. 
The  problem  at  this  point  breaks  down  into  two  sub-problems:  estimating  the  state 
using  Kalman-filter  theory  and  finding  the  control  signal  which  minimizes  the  cost  of 
equation  4.3.  Here  is  where  some  assumptions  are  made  to  simplify  the  design  and 
simulation  for  the  UWAA  planar  truss  reaction  mass  actuator  controllers. 

4.2.1  Optimal  Controller  Design 

For  the  application  of  the  optimal  control  design  techniques  to  the  UWAA  planar  truss 
RMA  controller,  several  assumptions  were  made.  First,  the  design  is  to  be  tested  only 
in  simulation  as  the  applicable  hardware  is  not  available  at  this  time  to  implement  a 
controller  with  its  own  dynamics.  Since  the  design  is  to  be  tested  only  in  simulation, 
it  was  assumed  that  all  the  states  were  available  for  feedback  thus  eliminating  the 
need  for  an  estimator  design.  In  practice,  this  could  be  accomplished  with  two  more 
accelerometers  mounted  at  the  truss  3/4  point  and  1/4  point  and  the  use  of  integrating 
circuits,  but  would  probably  prove  to  be  impractical.  However,  the  goal  of  the  research 
in  this  chapter  is  just  to  indicate  the  degree  of  performance  one  may  expect  from  the 
added  complexity  of  an  optimal  controller;  future  implementation  problems  will  need 
to  be  solved  later. 

The  second  assumption  is  to  ignore  noise  of  the  accelerometer  sensors  and  reaction 
mass  actuators.  Although  noise  problems  crept  up  in  the  experimental  application 
of  control  laws  in  chapter  three,  in  the  simulation  no  provisions  are  made  for  noise 
injection  to  the  sensors  and  actuators.  Thus,  we  may  measure  performance  gains  in 
eigenvalue  damping  but  lack  some  of  the  reality  noise  effects  bring  to  the  simulation. 
Without  the  inclusion  of  noise  the  LQG/Loop  Transfer  Recovery  technique  is  likewise 
not  available  as  a  design  tool.  Futtire  designs  for  an  optimal  controller  w’ill  have  to 
include  noise. 

Thus,  we  are  left  with  the  same  twelfth  order  plant  model  developed  in  sec¬ 
tion  2.3.1; 
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X  =  Ai2X  -f  B12U 

where  the  states  are  defined  as 


/  •  \ 

X\  ’ 

Htip 

X2 

Q3f4 

2:3 

Qmid 

X4 

9i/4 

^5 

^ta 

Qma 

X7 

Qiip 

93/4 

X9 

Qmid 

XlO 

91/4 

Xii 

9ta 

^  9ma  y 

The  control  inputs  are  likewise  given  to  be: 


u  = 


/  fb[tip) 
\  fb(mid) 


The  control  signal  will  be  given  by  the  feedback  law: 


u  =  —KcX 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


where  Kc  is  the  state-feedback  matrix  [13]  and  satisfies  the  deterministic  cost  func¬ 
tion  [12]: 


J  =  -  f  {z^Qz  +  iJ Ru)dt 

(4.9) 

2  Jo 

or 

1  roo  P  ^ 

(4.10) 

To  solve  for  K^.  it  is  recognized  that  equation  4.8  may  be  rewritten  as  [13]: 


u  -  -R-^lfP,x 


(4.11) 
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where  Kc  =  R  ^  Pc  and  Pc  satifies  the  algebraic  Ricacati  equation  [13] 

Pc  +  PcA  -  PcBR-^B'^ Pc  +  M^QM  =  0  (4.12) 

The  closed  loop  system  is  guaranteed  to  be  asymptotically  stable  for  a  stabilizable 
through  u  and  detectable  through  r  system.  Furthermore,  as  the  control  penalty 
weight  R  is  reduced  to  zero,  perfect  regulation  with  infinite  control  is  achieved  for 
our  minimum  phase  system  [12].  The  problem  of  design  rests  only  with  choosing  the 
proper  values  for  the  Q  and  R  weighting  matrices.  To  help  make  a  logical  choice  for 
the  weighting  matrices  Bryson’s  rules  were  applied  as  follows  [12]. 

Bryson’s  rules  for  selecting  the  diagonal  elements  Qa  and  Ra  are 


Q«  =  ^ 

(4.13) 

and 

Rii  = 

(4.14) 

imax 


where  Zimax  and  u^max  are  the  maximum  allowable  excursions  in  z,  and  [12].  The 
transformation  matrix  M  in  equation  4.4  was  chosen  to  be  equal  to  an  eigth  order 
identity  matrix,  thus  equating  z  with  the  states  x.  The  penalty  weights  are  thus 
associated  with  exactly  one  state  each  *.  The  weighting  values  were  chosen  by  choos¬ 
ing  the  maximum  excursions  for  truss  velocity  states  to  be  0.1  in/sec  and  to  be  0.01 
inches  for  the  truss  displacement  states.  The  reaction  mass  actuator  velocity  states 
were  chosen  to  have  limits  of  1.0  in/sec  and  limits  on  the  displacement  states  equal  to 


'  Another  scheme  was  successfully  chosen  to  pick  the  treinsformation  matrix  z  which  put  penalties 
on  the  truss  accelerations  at  the  four  truss  locations  and  on  the  two  RMA  displacement  states. 
The  transformation  matrix  had  the  form 


an 

012  •• 

01,11 

0l,12 

021 

022  •  • 

02,11 

02,12 

031 

032  •  • 

03,11 

03,12 

041 

042  •  • 

04,11 

04,12 

0 

0  • 

1 

0 

0 

0 

0 

1 

The  results  obtained  were  similar  to  those  obtained  using  the  states  as  the  penalized  variables  for 
gain  matrices  of  similar  magnitude. 
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the  maximum  throw  of  the  actuators  at  0.1  inches  of  travel  as  given  in  section  2.3.1. 
The  resulting  diagonal  elements  of  the  penalty  matrix  are  shown  in  table  4.1. 

Table  4.1;  Penalty  weights  on  the  states 


state 

Q„ 

^imax 

ESI 

state 

Qii 

K! 

9ttp 

Qn 

0.1 

100 

%ip 

Qn 

0.01 

10000 

93/4 

Q22 

0.1 

100 

^3/4 

Qsa 

0.01 

10000 

Qmid 

Q33 

0.1 

100 

0.01 

10000 

9i/4 

Q44 

0.1 

100 

9i/4 

<5 10,10 

0.01 

10000 

^ta 

Q55 

1.0 

1 

Qtip 

Qii.ii 

0.1 

100 

^tm 

Qee 

1.0 

1 

(Itip 

Qi2,\2 

0.1 

100 

The  penalties  on  the  control  inputs  were  chosen  initially  using  Bryson’s  rules  to 
limit  the  input  signal  to  the  amplifier  (set  at  unity  gain)  to  maximum  excursion  of  one. 
This  selection  of  a  penalty  resulted  in  an  abnormally  high  values  in  the  state  feedback 
gain  matrix  Kc-  After  some  trial  and  error  to  limit  the  magnitudes  of  the  individual 
feedback  gains  to  under  ten,  the  penalties  in  the  R  matrix  were  chosen  to  be  the  values 
presented  in  table  4.2  Solving  the  algebraic  Ricatti  equation  given  by  equation  4.12 

Table  4.2:  Penalty  weights  on  the  control  inputs 


input 

Rii 

^imax 

BUmI 

fb{tip) 

Ri\ 

0.0632 

250 

fb(mid) 

Rit 

0.0632 

250 

results  in  the  gain  matrix  /\c,  the  numerical  values  of  which  are  given  in  appendix  B.l. 
The  magnitude  of  the  largest  gain  is  less  than  ten,  which  when  the  accelerometer 
scale  factor  is  introduced,  is  under  100.  This  was  seen  as  a  reasonable  limit  for 
implementation.  The  problem  with  the  instrumentation  instability  of  section  3.5.2 
would  be  avoided  as  the  gains  assume  unity  gain  for  the  amplifier,  well  below  the 
instability  threshold  of  amplifier  gain  observed  at  gains  of  seven  or  greater. 
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4.3  Performance  of  the  Linear  Optimal  RMA  Controller  through  Simulation 

4.3.1  Linear  Frequency  Response 

The  RMA-optimal  controller  was  tested  using  the  linear  simulation  of  the  AJT  forced 
RMA-controlled  truss  model.  The  frequency  response  in  figure  4.1  shows  the  amount 
of  vibration  suppression  that  the  optimal  RMA  controller  introduced  to  the  system. 
Mode-1  suppression  was  8.67  dB,  mode-2  suppression  was  11.14  dB,  and  mode-3 
suppression  was  4.09  dB.  There  was  no  appreciable  suppression  of  the  the  fourth 
mode  vibrations  indicated  by  the  frequency  response  plot. 

4.3.2  Non-linear  Simulation  of  Modal  Time  Responses 

Modal  time  responses  of  the  truss  model  equipped  with  the  RMA  optimal  con¬ 
troller  were  generated  using  non-linear  simulation  of  the  system  by  driving  the  RMA- 
augmented  truss  model  with  the  AJTs  at  the  natural  frequencies.  The  optimal  con¬ 
troller  was  switched  on.  When  the  level  of  modal  excitation  was  approximately  the 
same  level  as  that  used  in  section  3.5,  the  AJTs  were  switched  off  and  the  linear 
simulation  decay  of  the  individual  modes  was  recorded.  The  eigenvalues  of  the  linear 
system  are  reported  in  table  4.3. 

Table  4.3:  Linear  model  modal  data:  optimal  RMA  controller 


frequency  (rad/sec) 

damping  ratio 

AC  vs.  open  loop 

mode  1 

9.19 

0.0193 

0.453 

mode  2 

58.387 

0.0837 

1.299 

mode  3 

145.56 

0.0158 

0.245 

mode  4 

260.27 

0.0187 

0.067 

RMA 

40.41 

0.10487 

0.146 

RMA 

43.19 

0.12170 

0.195 

Compared  to  the  classically  designed  RMA  controller  of  section  3.3.4,  the  damping 
of  the  first,  second,  and  third  modes  were  increased  by  a  greater  amount.  Mode-4 
damping,  however,  was  greater  for  the  classically  designed  controller.  The  optimal 
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frequency  (rad/sec) 


frequency  (rad/sec) 

Figure  4.1:  Linear  simulation  of  the  optimal- RM  A  controller  for  vibration  suppression 


controller  does  not  decrease  the  damping  in  any  of  the  modes  ^ls  was  the  case  with 
the  classical  controller.  However,  this  was  not  noticed  as  a  problem  for  the  classical 
controller  once  hybrid  actuation  was  introduced. 

4.4  l\on-Linear  Simulation  of  the  Hybrid  Classical- AJT/Optimal-RM A  Controller 


Hybrid  actuation  was  used  with  success  for  the  classically  designed  systems  of  chapter 
three  which  all  employed  velocity  feedback  without  any  crossfeeds.  The  ’piggy-back’ 
approach  running  the  AJT  controllers  cn  top  of  the  RMA  controlled  system  was 
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Figure  4.2:  Non-linear  simulation  modal  time  responses:  RMA  optimal  controller 


accepted  as  a  valid  approach  for  control.  The  question  remaining  was  if  hybrid  ’piggy 
back’  control  would  work  for  a  true  MIMO  system  with  crossfeed  terms. 

The  non-linear  simulation  was  performed  in  the  manner  common  to  this  thesis 
by  exciting  the  truss  model  with  the  AJTs  at  each  mode  with  the  RMA-optimal 
controllers  already  switched  on.  For  hybrid  control,  the  AJTs  were  switched  from 
exciting  the  truss  to  the  control  velocity  feedback  circuit  discussed  in  section  3.2.3 
while  the  RMA-optimal  controller  was  left  on  switched  on.  The  modal  time  responses 
are  presented  in  figure  4.3. 

The  mode-1  time  response  shows  a  significant  gain  over  the  independent  RMA- 
optimal  controller  system  in  .section  4.3.2.  The  three  remaining  modes  show  little 
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seconds 


Figure  4.3:  Modal  time  responses  hybrid 


seconds 
mode  4 


classical- AJT/optimal-RM A  controller 


change.  The  damping  ratios  were  calculated  and  are  presented  in  table  4.4.  The 
damping  ratios  do  indicate  that  the  hybrid  controller  only  significantly  effects  the 
first  mode.  The  mode-1  damping  ratio  is  0.1280,  well  above  the  elusive  10%  goal 
for  large  space  structures.  Although  the  simulation  generally  proved  to  have  better 
performance  than  the  actual  experimental  systems  in  chapter  three,  there  is  significant 
room  for  a  decrease  in  performance  should  the  optimal  controller  be  experimentally 
implemented. 
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Table  4.4:  Closed  loop  modal  data:  hybrid  classical- AJT/optimal-RM A  controller 


frequency  (rad /sec) 

damping  ratio 

AC  vs.  open  loop 

mode  1 

9.19 

0.1280 

1.275 

mode  2 

58.387 

0.0823 

1.292 

mode  3 

145.56 

0.0147 

0.213 

mode  4 

260.27 

0.0136 

-0.070 

4.5  Optimal  Control  Conclusions 

Optimal  control  design  for  the  RMAs  yields  a  system  which  does  not  sacrifice  mode- 
1  damping  when  damping  the  higher  modes  is  the  main  benefit  apparent  for  the 
optimal  Controller  of  this  chapter.  The  feedback  gains  in  this  section  are  larger  than 
the  gains  in  the  classical  system,  but  in  experimental  application  would  be  applied  in 
the  computer  hardware,  not  at  the  amplifier,  so  instrumentation/actuator  instability 
should  not  occur.  For  experimental  application,  an  estimator  design  and  the  inclusion 
of  sensor  and  actuator  noise  should  be  included  in  the  LQG  design  process. 

The  piggy-back  approach  is  again  a  valid  technique  for  application  of  hybrid  ac¬ 
tuation  using  velocity  feedback  AJT  controllers  on  board  the  optimally  designed 
RMA-controlled  planar  truss  for  vibration  suppression.  The  failure  of  the  hybrid 
controller  to  change  the  damping  of  the  second,  third,  and  fourth  modes  over  the 
optimally  controlled  system  is  disappointing.  Another  approach  might  be  applied 
in  lieu  of  the  piggy-back  approach  which  provides  for  non-linear  optimization  of  the 
AJTs  in  conjunction  with  the  RMAs. 


AJT 

non¬ 

linear 


approximate 

integrator 


Figure  4.4:  Non-linear  block  diagram  for  hybrid  classical- AJT/optimal-RMA  Con 
troller  (state  model  of  truss  in  section  A.2) 


Chapter  5 

CONCLUSIONS  AND  RECOMMENDATIONS 

5.1  Overview 

In  this  paper  the  set  up,  modeling,  and  control  of  the  UWAA  planar  truss  testbed  is 
described.  This  chapter  lists  the  conclusions  made  based  upon  the  work  performed  for 
this  study.  Recommendations  for  future  research  on  the  UWAA  planar  truss  testbed 
are  also  included. 

5.2  Conclusions 

The  purpose  of  this  paper  has  been  to  report  the  investigation  of  different  control 
laws  for  vibration  suppression  of  the  UWAA  planar  truss  via  the  available  means  of 
actuation,  the  air  jet  thrusters  and  the  reaction  mass  actuators.  The  following  tasks 
have  been  accomplished. 

•  Accurate  linear  and  non-linear  models  of  the  two  forms  of  space  realizable  actu¬ 
ation,  the  RMAs  and  the  AJTs,  have  been  formed  and  validated. 

•  The  existing  open  loop  truss  model  has  been  validated,  and  the  augmentation  of 
the  truss  model  with  the  RMAs  defined.  Comparison  with  the  USAF  Academy 
truss,  upon  which  the  UWAA  truss  is  based,  is  also  reported. 

•  An  array  of  classical  control  techniques  applied  to  the  vibration  suppression 
problem  was  investigated.  The  goal  of  10%  structural  damping  is  nearly  reached 
via  hybrid  control,  the  use  of  the  two  means  of  actuation  together. 

•  The  optimal  controller  design  problem  has  been  formulated.  Although  satisfac¬ 
tory  simulation  results  were  obtained  using  optimal-LQG  controlled  RMAs  in 
conjunction  with  classically  controlled  AJTs,  there  exists  the  need  to  account 
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for  the  non-linearities  in  the  optimization  problem.  The  optimal-control  solu¬ 
tion  of  this  paper  does  not  provide  for  much  higher  levels  of  damping  than  the 
classically  designed  hybrid  controller. 

•  The  set  up  of  the  working  truss  hardware  is  explained  and  the  modeling/hardware 
validation  wilt  serve  as  a  baseline  for  future  studies.  Computer  non-linear  sim¬ 
ulation  is  shown  to  accurately  reflect  the  results  achieved  using  the  hardware  in 
experimental  application.  The  computer  work  files  are  available  on  line  in  the 
truss  directory  the  University  of  Washington  [21]. 

5.3  Recommendations  For  Future  Research  on  the  UWAA  Planar  Truss 

The  next  step  for  future  study  should  be  the  implementation  of  the  control  laws 
reported  in  this  thesis  through  the  use  of  a  digital  controller  such  as  LABVIEW  or 
other  I/O  capable  software/hardware  packages  The  use  of  a  digital  controller  allows 
for  flexibility  in  implementing  several  different  control  laws  quickly  and  accurately, 
and  later  would  allow  for  the  implementation  of  high  order  optimal  linear  control 
laws  or  non-linear  control  laws. 

The  UWAA  truss  would  serve  as  an  ideal  platform  to  test  other  finite  element 
model  codes.  Modeling  of  the  non-linear  joints  of  truss  structures  along  with  the 
incorporation  of  both  viscous  and  coulomb  damping  terms  in  the  finite  element  code 
was  performed  by  W.  Keith  Belvin  of  Langley  Research  Center  [3].  The  joints  are 
reported  to  have  a  dead  band  (different  from  the  AJT  dead  band  of  this  thesis)  within 
which  a  joint  may  be  made  to  move  with  nearly  zero  applied  force.  This  modeling 
of  the  truss  indicated  that  joints  lower  the  overall  stiffness  of  the  truss  and  introduce 
additional  damping  terms,  substantially  affecting  the  low  frequency  vibrational  re¬ 
sponse.  A  non-linear  model  for  the  truss  would  be  an  ideal  intermediate  step  when 
progressing  from  the  testing  of  controllers  on  linear  models  to  actual  experimental 
implementation. 

A  parametric  study  testing  the  effectiveness  of  other  actuator  and  sensor  loca¬ 
tions  rather  than  the  truss  tip  and  midpoint  would  be  valuable.  Noting  that  it  is 
time  consuming  to  reconfigure  the  hardware,  most  of  the  study  could  be  performed 
via  simulation,  the  results  of  which  could  then  be  validated  through  the  use  of  the 
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hardware. 

Anothf :  RMA  might  be  built,  and  should  incorporate  such  features  as  user  vari¬ 
able  stiffness  and  damping  terms  through  the  use  of  electronic  feedback  instead  of 
the  present  preset  flexible  structure.  A  detailed  study  might  be  performed  where 
the  specifications  for  the  RMA  to  include  the  throw,  proof  mass,  range  of  required 
damping  and  stiffness,  and  location  (probably  tip  for  one  RMA)  are  all  identified 
before  the  construction  of  the  actuator  is  approved.  This  would  rectify  many  of  the 
limitations  of  the  current  RMAs. 

The  optimal  control  problem  for  the  RMA-augmented  truss  could  be  performed  in 
a  parametric  study.  LQ  techniques  which  include  the  state-estimator  design,  H2I Hiai 
methodology,  and  SANDY  design  all  need  to  be  explored.  Once  the  digital  controller 
is  operational,  all  of  these  controllers  would  be  easy  to  experimentally  implement  in 
the  truss  hardware. 

A  difficult  problem  would  be  to  develop  an  optimized  solution  for  the  non-linear 
AJTs.  An  easier  first  step  might  be  to  design  and  construct  an  AJT  with  linear 
characteristics  and  to  develop  the  optimal  control  problem  from  there.  True  optimal 
hybrid  control  would  thus  be  achieved. 

The  existing  truss  and  actuator  models,  hardware,  and  simulation  should  serve  as 
a  baseline  for  studies  in  advanced  hybrid  actuation  techniques  as  well  as  multi-input 
multi-output  control  law  studies.  The  truss  hardware  allows  for  the  incorporation  of 
active  truss  members  in  place  of  the  ordinary  structural  members.  Current  research  at 
NASA’s  Jet  Propulsion  Laboratory  in  Pasadena  [4]  is  being  conducted  utilizing  active 
structural  members  with  piezoelectric  actuators  and  colocated  force  and  vibrational- 
velocity  sensors.  The  control  laws  are  designed  to  control  the  mechanical  impedance 
of  the  truss  members.  The  mechanical  impedance  is  defined  as  the  Laplace  transform 
of  the  ratio  between  the  axial  force  applied  to  the  member  and  the  difference  between 
the  velocities  of  its  ends  and  is  analogous  to  electrical  output  impedance.  Ideally, 
creating  the  complex  conjugate  of  the  entire  structure’s  impedence  with  the  active 
members  through  the  use  of  control  laws  dissipates  the  energy  at  a  maximum  rate 
and  hence  provides  for  the  greatest  level  of  damping.  If  the  LIWAA  planar  truss  could 
be  made  to  shed  energy  in  this  way,  the  next  step  would  be  to  incorporate  hybrid 
control  by  including  the  AJTs  in  the  analysis.  Further  research  as  to  the  optimum 
placement  of  active  structural  members  could  also  be  performed. 
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Large  space  structures  will  have  an  additional  constraint  besides  lightweight  as  the 
cost  of  shuttle  flights  increase,  that  is  to  minimize  astronaut  extra- vehicular  activity 
(EVA)  time.  Research  ideas  which  involve  leading  edge  technologies  are  hypothesized 
here  to  meet  the  constraint  of  limited  EVA  time.  Trainable  alloys  which  change 
shape  when  a  small  cm  rent  is  applied  might  be  introduced  to  the  truss  structure. 
Lightweight  and  simple,  trainable  alloys  might  replace  the  current  structural  member 
actuators  used  today.  A  move  towards  embedded  sensors  in  the  structure  to  replace 
accelerometers  mounted  externally  would  again  simplify  construction  of  a  structure 
in  space  and  decrease  weight.  Lightweight  fiber  optic  relay  systems  would  increase 
reliability  of  sensor  signals. 

Finally,  an  intelligent  control  system  which  could  learn  characteristics  of  various 
truss  configurations  and  generate  controllers  would  replace  the  need  for  current  gain 
scheduling  routines  of  different  truss  configurations.  The  UWAA  planar  truss  hard¬ 
ware  is  fully  modular  so  that  appendages  could  easily  be  added  to  test  adaptive 
control  schemes. 
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Appendix  A 

TRUSS  LINEAR  STATE  SPACE  MODELS 


This  appendix  contains  the  numerical  values  for  the  two  truss  models  used  in  this 
thesis.  The  models  are  an  eighth  order  model  for  the  truss  without  actuators  attached 
and  a  twelfth  order  model  for  the  truss  with  RMAs  attached  to  the  truss  midpoint 
and  truss  tip. 


A.J  Open  Loop  Truss  Model  Without  Actuators  Attached 
The  open  loop  model  for  the  truss  dynamics  is  written  as 

i  =  Agx  +  Bgu  (A.l) 


(A.3) 
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The  outputs  were  chosen  beised  on  what  was  available  from  the  truss  hardware. 
For  the  accelerometers  mounted  at  the  tip  and  midpoint  measuring  the  lateral  ac¬ 
celerations  at  those  two  locations  the  output  state  space  equations  may  take  on  the 
form: 


y  = 


flu  012  ■  •  •  Ois 
031  <*32  •  •  •  038 


X  -t- 


0  0 
0  0 


u 


(A.4) 


The  numerical  values  which  comprise  the  A  and  B  state  matrices  are: 


As  = 


-1.4731 

3.2908 

-2.7590 

1.1269 

-11,358 

27,361 

-22,563 

9095.9 

1.1450 

-3.1484 

2.8300 

-1.2085 

9551.3 

-25, 162 

23, 264 

-9693.1 

-0.69798 

2.4674 

-3.5871 

2.4496 

-5689.0 

20,316 

-28,721 

20, 166 

0.18168 

-1.0198 

2.4646 

-3.3488 

1408.8 

-8189.2 

20,306 

-27,067 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

B8  = 


15.0207  0.75580 

-2.0949  -0.18716 

0.75580  6.3577 

-0.29383  0.04895 

0  0 

0  0 

0  0 

0  0 
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A. 2  Open  Loop  Truss  Model  with  RMAs  Attached  at  Midpoint  and  Tip 

The  open  loop  model  for  the  truss  dynamics  with  the  RMAs  attached  at  the  truss 
tip  and  truss  midpoint  is  written  as 


X  =  A12X  +  B12U 


(A.5) 


The  state  vector  is 


'  9Up  ^ 

X2 

?3/4 

xz 

Qmid 

X4 

Qi/4 

3^5 

^ta 

3^6 

Qma 

Xf 

Qtip 

Xg 

Qzli 

X9 

Qmid 

XlO 

9i/4 

Xii 

Qta 

K  ^12  > 

^  Qma  j 

(A.6) 


The  truss  stations’  velocities  and  displacements  have  the  same  variable  labels  as 
in  section  A.l.  The  two  new  variables,  and  qma  are  the  relative  displacements  at 
the  truss  tip  and  truss  midpoint  between  the  tip  and  midpoint  proof  mass  positions 
and  the  corresponding  truss  stations.  The  displacements  are  in  inches.  The  velocities 
are  in  inches  per  second. 

The  first  two  force  inputs  for  the  truss  tip  and  truss  midpoint  are  retained  since 
that  is  where  the  AJTs  were  located.  Inputs  three  and  four  are  the  forces  applied 
to  both  the  truss  station  and  the  proof  mass  at  each  of  the  respective  stations.  The 
input  vector  is  scaled  in  Ibf.  and  is  written  as 


'  ,  ' 

«2 

fmid 

«3 

fb{tip) 

^  fb{mid)  y 

i 

(A.7) 
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The  outputs  were  chosen  based  on  what  was  available  from  the  truss  hardware 
similar  to  what  is  reported  in  section  A.l.  For  full  state  feedback  an  identity  matrix 
was  used  for  C  and  a  zero  matrix  for  D  in  the  modeling. 

The  numerical  values  which  comprise  the  A  and  B  state  matrices  are: 

Ai2  = 

columns  1  through  6 


-1.4731 

3.2908 

-2.7590 

1.1269 

0.912440 

0.049127 

1.1450 

-3.1484 

2.8300 

-1.2085 

-0.12569 

-0.012165 

0.69798 

2.4674 

-3.5871 

2.4496 

0.045348 

0.41325 

0.18168 

-1.0198 

2.4646 

-3.3488 

-0.017630 

0.0031818 

1.4731 

-3.2908 

2.7590 

-1.1369 

-7.0818 

-0.049127 

6.9798 

-2.4674 

3.5871 

-2.4496 

-.045348 

-7.08032 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 
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columns  7  through  12 


-11,358 

27,361 

-22,563 

9095.9 

250.92 

13.604 

9551.3 

-25, 162 

23,264 

-9693.1 

-34.566 

-3.3688 

-5689.0 

20,316 

-28,721 

20,166 

12.471 

114.44 

1408.8 

-8189.2 

20,306 

-27,067 

-4.8482 

0.88110 

11,358 

-27,361 

22,563 

-9095.9 

-1947.5 

-13.604 

5689.0 

-20,316 

28,721 

-20, 166 
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-1961.5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

15.0207 

0.75580 

-15.207 

-0.75580  ' 

-2.0949 

-0.18716 

2.0949 

0.18716 

0.75580 

6.3577 

-0.75580 

-6.3577 

-0.29383 

0.04895 

0.29383 

-0.048950 

-1.5207 
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Appendix  B 

FULL  STATE  FEEDBACK  GAIN 

B.l  OptimaJ  Control  Gaia  Kc 

1.5523e  -  1  2.0046e  -  1 

6.9095e-l  1.8848e-l 

4.1918e-l  -1.8663e-l 
3.6230e-2  -2.9339e  -  2 
1.0046c -  1  -4.4728c -3 
-3.9297c  -  3  5.5639c  -  2 

-4.6436c  00  3.6697e  00 

1.7244c  00  -5.1648e  -  1 
7.7796e-2  -8.6408e  00 
7.4760c -2  -2.6993e  00 
1.2032e-2  -5.3007c  -  2 
5.7819c -2  1.1015c -2 


Appendix  C 

EQUIPMENT  LIST 

This  appendix  contains  a  partial  listing  of  the  equipment  used  in  the  research 
reported  in  this  thesis.  The  list  is  by  no  means  all-inclusive  but  gives  a  good  idea 
as  to  what  hardware  is  required  to  set  up  the  planar  truss  experiment.  Manufactuer 
information  are  included  in  boldface  for  applicable  items. 

C.l  Truss  Hardware 

•  Meroform  522117  tube-m  0.353  22M12  alum  E6-CO  (61) 

MERO  Structures,  Inc. 

N112  W1880  Mequon  Road 
Germantown,  WI  53022 
phone:  (414)255-6553 

•  Meroform  522218  tube  0.500-m  22M12  alum  E6-CO  (20) 

•  Meroform  510018  std.  node  M12  Chrome  (42) 

•  3/8  in-dia  steel  ball  bearing's  (40-1-6  for  RMAs) 

•  Level  table  4’x8’  top  (3) 

C.2  Custom  Fabricated  Truss  Hardware 

•  Steel  table  for  cantilever  support 

•  Cantilever  mounting  point  (box  beam) 

•  Steel  battens  (20) 
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•  1  7/8  X  4  in  stainless  steel  type  301  sheets  (40+) 

•  Accelerometer  housings 

C.3  Air  Jet  Thruster  Hard, /are 

•  Solenoid  valve  MAC;113B551BAAA  (4) 

Fiero  Fluid  Power  Inc. 

10515  East  40th  Ave. 

Denver,  CO  80239 
phone;  (303)  373  -  2600 

•  Polyvinyl  tubing  1/4  in-inner  dia  3/8  in-outer  dia  (80  feet) 

•  Brass  nozzle  (hose  barb)  1/4  in  i.d.  x  1/4  in  in.  MIP  (12) 

•  65-psig  regulated  compressed  air  supply 

C.4  Custom  Fabricated  Air  Jet  Thruster  Hardware 

•  Manifold  (for  delivery  of  compressed  air  to  individual  AJTs) 

C.5  Reaction  Mass  Actuator  Hardware 

•  LDS  V102  vibration  generator  (2) 

Ling  Dynamic  Systems,  Inc. 

60  Church  Street 
Yalesville,  CT  06492 
phone:  (203)  265  -  7966 


•  Swivel  linkage  (stinger)  (2) 
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C.6  Custom  Fabricated  Reaction  Mass  Actuator  Hardware 

•  Aluminum  calibration/test  jig 

•  Aluminum  proof  mass 

•  Aluminum  cradle  and  grooved  track 

•  Numerous  adaptors  for  motor /stinger /node/ test  jig  attachments 

C.7  Instrumentation  and  Electronics 

•  Comdyna  GP-6  Analog  Computer  (2) 

Comdyna,  Inc. 

305  Devonshire  Road,  Barrington 
IL  60010 

phone;  (708)  381  -  7560 

•  Hewlett  Packard  HP35665A  Dynamic  Signal  Analyzer 

•  MDB-5  load  cell  and  TM-2  amplifier 

Transducer  Techniques 
43178  Bussiness  Park  Drive  B-101 
Rancho  Temecula,  CA  92390 
phone:(714)676-1200 

•  Sundstrand  QA-700  accelerometer  (4) 

•  Control  circuit  hardware  (misc.  see  also  Schulthess  [17]) 

•  PA  11 503  ELECTRO- MIKE  position  sensor 

ELECTRO  Corp. 

1845  57th  Street 
Sar2isota,  FL  33580 


